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A  "split-field"  infrared  optical  system  has  been  designed  and  constructed  to 
provide  simultaneous  image  pairs  in  a  single  frame  of  an  infrared  (IR)  imager, 
differing  only  in  the  direction  of  linear  polarization.  The  optical  train  is  afocal, 
allowing  its  use  with  a  variety  of  infrared  imaging  devices.  The  system  can  operate 
in  both  long-wave  IR  (  8-12  um)  and  mid-wave  IR  (3-5  urn)  with  interchangeable 
polarizing  splitter  plates.  Previous  work  at  the  Naval  Academic  Center  for  Infrared 
Technology  (NACIT)  at  the  Naval  Postgraduate  School  has  demonstrated  that 
significant  image  improvement  in  infrared  image  contrast  can  be  obtained  by  use 
of  polarization  filtering,  especially  for  targets  at  sea  through  suppression  of 
polarized  sea  background  emission  in  the  long  wave,  or  of  polarized  reflection  in 
the  midwave.  That  work  utilized  digital  subtraction  of  sequential  image  pairs  with 
orthogonal  polarizations,  but  suffered  from  inability  to  obtain  simultaneous 
images  and  also  from  problems  due  to  reflections  from  external  polarization  filters. 
Both  of  these  problems  are  eliminated  with  the  new  split-field  technique.  Prelim- 
inary tests  of  this  system  with  an  AGA-780  imager  were  carried  out  in  a  field 
experiment  using  an  Arleigh  Burke  DDG  class  ship  as  target  at  varying  distances 
at  sea.  For  comparison,  images  were  also  taken  in  successive  pairs  in  time  with  the 
AGA-780  viewing  the  scene  directly  (  without  the  split-field  adjunct)  with 
interchangeable  internal  polarizing  filters  in  a  rotatable  filter  wheel.  Subjective 
and  numerical  analysis  of  the  data  from  the  field  experiment  demonstrated  good 
image  quality  and  contrast  improvement,  and  potential  for  future  utility. 
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I.   INTRODUCTION 

Previous  work  done  at  the  Naval  Postgraduate  School  has 
shown  that  infrared  image  contrast  for  sea  targets  can  be 
improved  by  polarization  filtering  (Refs:  20,21,22). 
Unwanted  background  radiance  which  clutters  the  image  can 
consist  of  emission  and  reflection  components.  At  grazing 
angles  an  appreciable  part  of  the  water  surface  emission  is 
linearly  polarized  with  the  electric  vector  vertical,  while 
direct  or  indirect  solar  radiation  reflected  from  the 
surface  shows  horizontal  polarization.  Other  elements  of  the 
scene,  such  as  the  sky  and  man-made  objects,  often  emit  or 
reflect  radiation  that  is  unpolarized  or  polarized  in  a 
different  direction  from  that  of  the  sea  surface.  As  a 
result,  viewing  the  scene  through  a  polarizing  filter,  or 
subtracting  two  otherwise  identical  images  except  for 
orthogonal  polarization  directions,  can  improve  the  target 
contrast  against  background. 

Polarization  filters  for  the  infrared  are  only  recently 
available.   These  consist  of  "wire-grid"  structures 
deposited  by  photolithography  and  etching  on  transparent 
substrates.  Such  substrates  also  need  to  be  antiref lection 
coated  to  reduce  unpolarized  reflection  by  the  dielectric 
substrate. 


Special  problems  arise  when  dealing  with  IR  as  opposed 
to  visible  radiation.   Because  all  materials  emit  radiation 
at  temperatures  above  0  degrees  K,  sensor  background 
radiation  often  contaminates  the  final  image.   Very  low 
temperatures  must  be  used  at  all  points  from  which  radiation 
can  enter  the  detector  in  order  to  reduce  the  amount  of 
unwanted  radiation  entering  below  acceptable  levels.   Liquid 
Nitrogen  and  thermo-electrically  cooled  bodies  are  commonly 
utilized.   Also,  the  'narcissus  effect'  can  occur  in  the 
final  image  when  the  cold  detector  'sees'  a  reflection  of 
itself  from  a  flat  or  back- focussing  optical  surface  (Ref. 
23) .  Various  methods  to  eliminate  these  problems  are 
implemented  in  the  techniques  presented.   These  are  key 
driving  factors  in  the  utilization  of  internal  filtering  and 
in  the  design  of  the  Naval  Postgraduate  School  Split-Field 
infrared  polarized  imaging  system. 

Important  to  the  study  of  image  improvement  is  the 
quantitative  characterization  of  image  intensity 
differences.   In  previous  experiments  done  at  the  Naval 
Postgraduate  School  in  which  the  orientation  of  the 
polarizer  had  to  be  changed  manually  between  recording 
orthogonally  polarized  images,  there  were  changes  in  the 
image  in  the  time  interval  during  this  process 


(Ref s . 1, 21, 22) .  Consequently  it  was  impossible  to  have  two 
identical  images  for  comparison.   To  surmount  this  problem  a 
Split-Image  optical  system  was  designed  and  built  at  NPS.  In 
this  system  two  orthogonally  polarized  images  are  recorded 
simultaneously.  The  images  can  then  be  digitized  and 
mathematically  subtracted  to  demonstrate  the  contrast 
improvement.   Specific  design,  implementation,  performance 
and  image  analysis  will  be  explained  and  presented. 

A  description  is  also  given  of  adaptation  of  the  AGA 
imager  with  polarizing  filters  installed  in  the  internal 
filter  wheel.  This  system  is  useful  in  recording  sequential 
images  with  orthogonal  polarizations. 

This  thesis  explores  the  phenomenon  of  infrared 
polarized  radiation,  background  and  target  polarization 
characteristics,  and  experimental  techniques  for  filtering 
and  improving  image  quality. 


II.  BASIC  INFRA-RED  THEORY 

To  fully  understand  the  phenomena  associated  with 
imaging  of  polarized  infrared  radiation  it  is  necessary  to 
understand  the  basic  physics  used  to  characterize  infrared 
phenomena.   This  chapter  deals  with  the  laws  of  thermal 
radiation,  radiation  measurement  metrics,  and  discussion  of 
thermal  imaging  fundamentals.   These  short  presentations 
attempt  to  give  some  insight  into  the  selection  of  the 
infrared  band  and  how  it  is  used. 
A.    INFRA-RED  TERMINOLOGY  AND  UNITS 

In  the  study  of  radiation,  the   radiant  power  exchanged 
between  the  target,  medium,  and  detector  is  of  interest.   A 
set  of  parameters  is  defined  to  describe  the  radiation  at 
each  point  of  interest.   Table  2.1  shows  some  of  the 
parameters  and  their  units  of  measurement  (Ref .  19) .   The 
last  column  in  Table  2.1  denotes  the  units  in  the  metric 
system.   The  entire  content  of  this  system  of  units  will  not 
be  addressed,  but  some  particular  units  need  to  be 
emphasized.   Temperature  is  measured  in  Kelvins  (K) .   The 
steradian  (Sr)  unit  of   solid  angle  is  derived  from  defining 

the  total  solid  angle  around  the  center  of  a  sphere  as  47i 

steradians.   Note  that  these  quantities  are  commonly  used  in 
inconsistent  units.   This  is  done  to  match  imaging  system 
dimensions  more  closely  for  system  performance  calculations. 


Symbol 

Name 

Description 

Units 
(metric) 

M 

Radiant 

Radiant  Flux  leaving 

Watts/cm2 

Flux 

an  element  of  a 

Density 

surface  per  unit 

(Exitance) 

area. 

E 

Irradiance 

Radiant  power  per 
unit  area  incident  on 
a  surface. 

Watts/cm2 

I 

Radiant 

Radiant  power  leaving 

Watts/Sr 

Intensity 

a  point  source  per 
unit  solid  angle. 

L 

Radiance 

Radiant  power  leaving 
or  arriving  at  a 
surface  at  a  point  in 
a  given  direction  per 
solid  angle  per  area 
normal  to  that 
direction. 

Watts/cm2Sr 

Table  2.1   Radiometric  Units 


All  other  physical  parameters  are  measured  in  conventional 
units.   Along  the  imaging  path,  M   characterizes  the  power 
leaving  the  target  surface  per  unit  area.   Some  of  this 
radiation  is  emitted  by  the  target  and  some  is  reflected  due 
to  Irradiance   from  another  source.   The  distant  transport  of 
the  radiant  energy  from  the  target  is  described  by  the 
Radiant   Intensity,    I,  which  includes  the  geometrical 
spreading  with  range.  I   is  the  most  useful  parameter  in 
calculating  the  maximum  range  of  an  imaging  detector.  When 


the  radiation  reaches  the  detector,  the  Irradiance   at  the 
receiver  of  the  target  radiation  is  converted  into  a  signal 
for  recording. 
B.    LAWS  OF  THERMAL  RADIATION 

In  order  to  detect  a  target,  it  is  necessary  to 
characterize  the  radiation  coming  from  the  surface.   Most 
infra-red  detectors  only  operate  efficiently  in  a  small  band 
of  wavelengths,  so  the  thermal  radiation  from  a  target  must 
be  known  in  that  range.   A  set  of  well  known  radiation  laws 
exists  that  describe  the  spectral  emittance  of  thermal 
sources . 

1 .    Black  Body 

A  black  body  is  known  simply  as  an  ideal  emitter  and 
absorber.   By  Kirchoff's  Law,  for  a  black-body  the  spectral 

absorptivity (a)  =  emissivity (e)  =  1,  so  the  reflectance (p) 
and  transmissivity  (x)  must  be  =  0 .   For  a  black  body  this 

is  true  at  all  wavelengths.  The  spectral  radiant  emittance 
curve  for  a  black  body  is  shown  by  Fig.  1,  and  is  described 
by  Planck's  Law. 


25 


30 


Figure    1 
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2 .   Planck  f  s  Law 

Planck's  Law  is  the  distribution  shown  in  Fig.  1.   It 
is  given  by 


M^iT) 


_  2nc2h/K5 


he 
KkT 


(e  ^-1) 


{N/  cm2]im) 


where  h  is  Planck's  constant,  and  k  is  Boltzman's  constant. 
The  spectral  radiance  of  the  black  body  can  be  reproduced  by 
the  Planck  distribution.   For  non-black  body  emitters,  the 


radiant  exitance  is  reduced  by  the  emissivity (s)  .  The 


radiant  exitance  is  then 


^=sx^ 


lack-bodv 


2) 


with  8<1  for  all  non-black  bodies.  Examples  of  emissivities 


Material 

Emissivity 

Polished  Al 

.05       (373K) 

Steel  Oxidized 

.79       (573K) 

Water (frost/ice) 

.98/. 96   (263K) 

Table  2.2   Emissivities  of  common  materials 

of  practical  materials  are  shown  in  Table  2.2  (Ref.  19) 

If  6X(T)  =  constant  independent  of  wavelength,  then  the  body 

is  called  a  grey  body.   For  most  surfaces  this  is  not  the 
case,  and  the  spectral  emissivity  varies  greatly  with 
wavelength.   Many  materials  behave  as  grey  bodies  over  a 
limited  range  of  wavelength. 

3.   Wien  Displacement  Law 

In  Fig.  1,  a  maximum  radiant  exitance  occurs  at  a 
wavelength  A^^.   A^ax  can  be  found  by  differentiating  the 


Planck  distribution  function  and  equating  to  zero.   The 


result   is   the  Wien   Displacement   Law. 

X     T=2898. 18]imK  (3) 

max  ~ 

Using  Wien's  Displacement  Law,  the  wavelength  with  the 
highest  spectral  radiance  is  known  with  the  surface 
temperature . 

4 .   Stefan-Boltzman  Law 

If  the  radiant  exitance  is  integrated  over  all 
wavelengths  for  a  given  temperature,  the  total  power  emitted 
per  unit  area  can  be  found.   The  equation  is  found  to  be, 
for  a  black  body 

M{T)=OTi  (4) 

ct  is  the  Stefan-Boltzman  constant,  a=5.67  x  10"12  W/cm2K4. 

For  other  bodies,  application  of  the  spectral  emissivity 

must  be  done  over  all  wavelengths. 

C.   ELECTROMAGNETIC  BANDS  AND  ATMOSPHERIC  WINDOWS 

To  understand  why  imaging  is  performed  in  the  infra-red 
band,  knowledge  of  the  bands  and  the  atmospheric  attenuation 
'windows'  is  necessary.   All  bodies  radiate  thermally  at  all 
wavelengths,  so  it  would  appear  logical  to  detect  the 
radiation  at  the  peak  wavelength  defined  by  Wien's 
Displacement  Law.   This  is  not  done  in  every  case,  due 
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partially  to  atmospheric  effects  on  the  emitted  radiation. 
A  brief  discussion  of  how  the  target  and  atmosphere  interact 
follows . 

1 .   Electromagnetic  Radiation  Bands 
The  infrared  region  of  the  electromagnetic  spectrum  lies 

between  the  radio  (longer  wavelength)  and  the  visible 

(shorter  wavelength) regions .   Fig.  2  shows  this  region  of 

the  spectrum  (Ref.  19).   Its  bandwidth  is  approximately 

.l\m   to  lOOOjim.   The  longest  wavelength  of  interest  for 

conventional  IR  imaging  is  approximately  15(am  since  most 
targets  and  modern  quantum  detectors  operate  below  this 
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Figure  2.   Partial  Electromagnetic  Spectrum 
wavelength.   termed  the  far  infrared  limit.   The  region 
nearest  the  visible  is  the  near  infrared.   The  "millimeter 
wave"  band  occurs  around  lOOOjxm. 
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2.   Atmospheric  "Windows' 

The  processes  that  effect  the  radiation  on  its  path 
through  the  atmosphere  are  extremely  complex  and  will  only 
be  considered  briefly.    Attenuation  affects  the  amount  of 
radiant  power  incident  on  the  optical  system,  which  in  turn 
affects  the  overall  performance  of  the  system.   The 
attenuation  mechanisms  of  primary  concern  are  absorption  and 
scattering  of  the  radiation.    Limiting  or  avoiding  the 
amount  of  absorption  and  scattering  in  the  target  path  is 
essential  for  any  thermal  imaging  system  because  of  the 
severity  of  the  effects.   Suspended  aerosol  particles  are 
the  primary  contributors  to  scattering  of  the  radiation  in 
the  IR.  Aerosols  include  dust,  fog,  salt,  ash,  and  many 
other  airborne  particles.   Molecules  such  as  H20  and  CO,  are 
the  predominant  absorbing  species  in  the  atmosphere  (Ref. 
19)  and  define  the  windows  in  which  IR  radiation  can  pass 
through  the  atmosphere  with  any  effectiveness.   This  is  due 
to  the  structure  of  the  molecules  and  the  modes  in  which 
they  can  be  excited  and  thus  absorb  energy.   Figure  3  shows 
a  plot  of  atmospheric  transmittance  versus  wavelength  over 
the  near  to  far  infra-red.   This  plot  was  generated  by  an 
atmospheric  propagation  modelling  code,  LOWTRAN  7.   LOWTRAN 
7  has  resident  atmospheric  models  for  aerosol  and  humidity 
content  in  addition  many  other  atmospheric  characteristics. 
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The  model  used  for  the  plot  is  for  a  Mid-Latitude  Summer 
atmosphere,  5km  path  at  sea-level.   The  'Rural'  designation 
is  for  the  aerosol  model  describing  an  air  mass  over  land 
with  primary  components  being  dust,  sand,  and  hygroscopic 
particles  [19]  . 
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Figure  3.   Atmospheric  Windows 

Transmission  with  no  attenuation  from  the  figure  has  a 
transmittance  value  of  1.   The  windows  are  depicted  by  the 
regions  of  maximum  transmittance.   It  is  clear  that  any 
imaging  would  need  to  take  place  for  radiation  found  in 
these  regions .   Commonly  used  wavelengths  for  imaging 

military  targets  are  around  4 . Oum  (short  wave)  for  hot 

sources  such  as  exhaust  plumes  and  power  plant  components, 

and  8-12|im  (long  wave)  for  surfaces  nearer  to  ambient 

temperatures.   The  relation  of  the  temperature  to  desired 
wavelength  band  follows  the  Planck  distribution  for  maximum 
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radiated  power.   This  correlation  between  radiated  power 
wavelengths  and  propagation  windows  must  be  considered  when 
using  IR  radiation  for  imaging  is  desired. 
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III.   POLARIZATION  PHYSICS 

Polarization  occurs  in  the  reflected  and  emitted  IR 
radiation  of  targets  and  backgrounds.   For  purposes  of 
enhancing  target  contrast  (as  exemplified  by  ship  targets  at 
sea) ,  it  is  necessary  to  determine  the  degree  of 
polarization  present  and  its  orientation.   Fundamental 
physical  laws  and  polarization  equations  describe  the 
important  parameters  of  polarization  of  the  target,  and  sea 
and  sky  background.   Quantification  of  these  parameters 
leads  to  expected  increases  in  contrast  and  is  necessary  in 
designing  and  implementing  polarizing  imaging  systems. 
A.    STOKES  PARAMETERS 

Sir  George  Gabriel  Stokes  (1852)  discovered  that  the 
polarization  state  of  any  light  source  could  be  completely 
described  by  four  observable  quantities  which  have  become 
known  as  the  Stokes  Parameters.   This  is  done  by  taking  the 
time  average  of  the  E  field  vector  amplitude  described  by 
the  polarization  ellipse,  thus  yielding  an  observable 
intensity  (Ref.17).   The  resultant  four  parameters  are 


I=E2+E2 

s  p 

m=e2-e2  (5: 

s        p 

C=2E  E  cos5 

s     p 

S=2E  E  sin5 

s     p 


Es   and  Ep   are   the   normal    and  parallel   components   of   the  wave 
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electric  field,  with  respect  to  the  plane  of  incidence. 
For  a  horizontasl  surface  the  vertical  component  (p)  lies  in 
the  plane  of  incidence,  also  referred  to  as  parallel 
orientation.   J  denotes  the  total  intensity  of  the  light,  M, 
the  intensity  of  linear  fs'  or  'p'  polarized  light,  C,  the 
linear  ±45  degree  polarized  light,  and  S,    the  circularly 
polarized  intensity.   8  is  the  phase  difference  between  the 

s  and  p  components  of  the  E  field.   The  degree  of 
polarization  is  the  ratio  of  the  intensities  of  the 
polarized  to  the  total  intensity  [3] . 


1/2 

(M2+C2+S2) 
P= 
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The  polarization  of  the  radiation  to  be  filtered  comes  from 
effects  due  to  emission  and  reflection.   Specific  laws  have 
been  developed  to  describe  these  phenomena  and  their 
contribution  to  the  Stokes  parameters  will  be  shown. 
B.    EMISSION  POLARIZATION 

The  target  and  background  radiance  produced  by  thermal 
emission  from  surfaces  is  inherently  polarized  [9].   Of 
specific  interest  for  imaging  targets  at  sea  is  the  emission 
polarization  of  the  sea  and  sky  background  since  the  primary 
contrast  improvement  will  be  obtained  by  filtering  of  the 
background. 
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1.    Sea  Emission  Polarization 

To  find  the  emission  polarization  of  the  sea,  the 
surface  can  be  treated  as  an  interface  with  coordinates 
shown  in  Figure  4  (Ref .  21) . 


Sun 


Wir.rl 


Figure  4.   Coordinate  System  of  Sea  Surface 


Y  is  the  azimuthal  direction  of  the  sun.   S  and  D  are  the 
directions  of  the  sun  and  observer,  respectively.   N  is  the 

wave  facet  normal,  0  and fy   the  zenith  and  azimuth  angles  for 
each  vector.   The  light  coming  from  the  upper  side  of  the 
interface  (+z)  consists  of  a  reflected  and  an  emitted 
component,  such  that  the  total  radiance  is  given  by 
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N(Q)=Np(Q)+N£  (9) 


The  light  emitted  from  the  surface  can  be  considered  as 
radiation  coming  from  the  sea  water  side  of  the  interface 
undergoing  refraction  at  the  surface.   Behavior  at  the 
interface  can  be  derived  from  the  Fresnel  equations  and  is 
discussed  in  this  fashion  by  Sandus  (Ref.  9)  and  Sidran 
(Ref.  11)  following  Reference  16.   The  emissivity  of  the 
surface  is  given  by  Reference  16. 


£  =  ■ 


(l+O 


zr  is  dependent  on  the  orientation  of  the  polarization 
(parallel  or  perpendicular)  where: 
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s+cos20 
z±= 

y2cos6[  (s+a-sin26)1/2] 

s+(a2+jb2)cos26 
z- 

y/2cose[s(a2+b2)  +  (a-sin2®)  (a  2-b2)  +  (2ab2)  ] 1/2 

a=Re{n2) 
b=Im(n2) 

s=J  (a-sin26)2+jb2 

n  is  the  complex  index  of  refraction  of  the  medium,  0,  the 

angle  of  emission.   These  equations  hold  true  for  any 
specular  surface.   Dependence  on  wavelength  comes  from  the 
index  of  refraction  (Ref.  11).   This  analysis  indicates  that 
both  the  emitted  and  reflected  components  of  surface 
radiance  will  show  polarization,  predominately  in  orthogonal 
directions.   Note  that  consideration  is  given  only  to  the 
linearly  polarized  light.   For  most  materials,  there  is  only 
a  small  range  of  angles  at  which  elliptical  polarization 
occurs.   Even  in  these  cases  the  phase  difference  between 
orthogonal  components  is  small  enough  (due  to  small 
dielectric  constants)  that  an  approximation  of  linearly 
polarized  light  at  all  angles  is  valid  (Ref.  3) .   This 
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theory  has  been  applied  to  the  sea  surface  and  proven  by 
Sidran  (Ref.  11),  Baesner  and  McCoyd  (Ref.  8),  and  Gregoris, 
Yu,  Cooper,  and  Milne  (Ref.  20).. 

In  dealing  with  the  sea  surface,  wave  structure  alters 
the  coordinates  of  interest  in  Fig.  4.   The  wave  structure 
is  considered  to  compose  a  distribution  of  emitting  facets 
with  unit  direction  toward  the  observer  defined  by  the 
vector  s     (see  Fig.  4),  where 


s-s1i+s2j+s3£ 


10 


The  orientation  of  the  normals  to  the  facets  with  respect  to 
the  normal  to  the  plane  defined  by  §    and  the  z  axis  is 
defined  by  the  angle  a.   The  effective  emissivities  of  each 

facet  referred  to  the  vertical  and  horizontal  directions, 
are  given  by 


e'  =e  cos2a+e,sin2a 

(11 


e'=e  sin2a+8,cos2a 


The  prime  symbol  denotes  the  emissivity  in  the  observer 
direction  relative  to  the  vertical  to  the  plane  containing 
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the  z  axis  and  the  observer  direction  §     (Ref .  8) .   The 
distribution  of  the  slopes  of  the  facets  has  been 
experimentally  determined  by  Cox  and  Munk  (Ref.  15) .   The 
Gram-Charlier  equation  with  empirically  derived  constants 
adequately  describes  the  distribution  of  wave  slopes  (Ref. 
15) .   The  resulting  wave  slope  distributions  are  considered 

valid  for  wind  speed  <14m/sec  and  wave  slope  angle  (P)<23 
degrees.   Plots  of  the  Cox  and  Munk  distribution  are  shown 
in  Fig.  5  from  (Ref.  11).   Because  shadowing  occurs  on  the 
back  surfaces  of  the  waves  in  the  direction  of  the  observer, 
a  correction  should  be  applied  to  the  equations  as  used  by 
Basener  and  McCoyd  (Ref.  8)  for  percent  polarization 
calculation.   Saunders  (Ref.  10)  gives  a  correction  to  the 
Basener  and  McCoyd  results.   Figure  6  shows  the  polarization 

of  emission  for  1 .  Ojxm  wavelength  and  14m/s  wind  speed  after 

Sidran  (Ref.  11).     The  percent  of  polarization,  Q  is 
defined  by 

c=ioo.i!^  (12) 

± 


The  results  in  Fig.  6  hold  true  for  all  wavelengths  <  lOum 
excluding  the  absorption  bands  (Ref.  11) .   Note  that  the 
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largest  polarization  effects  are  seen  at  grazing  incidence. 
2.   Sky  Emission  Polarization 

The  IR  radiation  from  the  sky  is  unpolarized  (Ref .  2) . 
Use  of  polarized  filtering  components  will  not  contribute  to 
the  contrast  improvement  measurements  when  the  sky  is  viewed 
directly.   The  sky  emission  in  the  image  will  be  subject  to 
attenuation  by  the  polarizing  filter  to  the  degree  of  the 
unpolarized  transmissivity  of  the  filter.   Fig.  7  (Ref.  10) 
shows  the  sky  radiance  as  a  function  of  zenith  angle. 
Specific  radiance  characteristics  can  be  seen  in  Fig.  7. 
When  viewed  toward  the  zenith,  the  sky  emits  less  due  to  the 
fact  that  the  atmosphere  consists  of  water  vapor  and 
aerosols  (absorbers)  which  also  act  as  good  emitters,  but  is 
less  dense. 
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Figure  5.   Facet  slope  angle  distribution 
(a) crosswind,  (b)  upwind-downwind  [11]. 
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Figure  7.   Sky  Radiance  as  a 
function  of  zenith  angle  from 
Saunders  [10] . 
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Sky  radiance  increases  with  zenith  angle  toward  the  horizon 
to  that  of  a  black-body  with  atmospheric  temperature.   This 
is  due  to  the  maximum  thickness  of  the  atmosphere 
encountered  at  that  angle  and  even  the  small  emissivities 
of  the  absorbers  combine  (Ref .  10) . 

2.    Target  Emission  Polarization 

Although  primary  contrast  improvement  in  the  at-sea 
scenario  will  be  gained  by  filtering  of  sea  emission 
polarization,  there  is  polarization  due  to  emission  from  the 
target  itself  that  must  be  considered.   Painted  surfaces 
such  as  are  found  on  all  types  of  vehicles  and  ships  display 
emission  polarization  (Ref.  3,4).   Figure  8  shows  21 
measurements  taken  on  various  paint  samples  ranging  from 
very  smooth  torough  sand-paint  mixtures  (Ref.  4) . 


—  Ideol  Surfoce 
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Measured  polarized  signature  components  in  the  7.5-12  micron  band 
for  selected  paint  samples  viewed  at  45°  from  normal. 


Figure    8.    Ref    [4] 
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The  line  for  an  ideal  surface  in  Fig. 8  denotes  the  amount  of 
polarization  that  would  be  seen  from  a  perfectly  specular 
reflector.   The  rougher  samples  show  less  polarization  than 
the  smoother.  The  amount  of  polarization  depends  on  the 
index  of  refraction  of  the  paint  and  any  degradation  due  to 
surface  roughness. 

C.    REFLECTION  POLARIZATION 

For  thermal  radiation  at  a  surface 

p(0)+a(e)+T(6)=l     e(e)=a(0)         (13) 

As  stated  previously,  the  total  amount  of  energy  at  a 
surface  consists  completely  of  the  reflected,  absorbed,  and 
transmitted  fractions.   For  a  perfectly  opaque  body, 
therefore 

e(6)+p(6)=l  (14) 

The  emissivity  (8)  was  treated  in  the  previous  section;  the 

reflectance  (p)  of  the  sky,  sea  surface,  and  target  must  now 
be  considered.   Absorption  will  have  an  effect  on  the 
reflectance  for  each  surface  considered  and  will  act  to 
degrade  the  reflectance. 
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1 .    Sea  Reflection  Polarization 

The  application  of  the  Fresnel  Equations  to  a  Cox-Munk 
sea  surface  slope  distribution  can  be  shown  to  predict 
polarization  in  reflectance  as  well  as  emittance  with  the 
following  relations [11] . 


e±-i-Pi   e-i-P| 


e=  1/2  (£±+  £ 


p=l/2(px+p| 


(15) 


£  =  1-p 


The  sea  surface  acts  as  an  opaque  body.   Values  of  s  may  be 

obtained  as  discussed  in  the  previous  section.   The 
percentage  polarization  seen  due  to  reflection  is  defined  as 

[11] 


P  "Pn 

P=100-^^L 


P.+  P. 


2 .    Sky  Reflectance 

As  stated  in  the  emission  section,  the  sky  acts  as  a 
poor  emitter  near  the  zenith  angle  due  to  the  low  optical 
depth.   Because  of  the  high  content  of  atmospheric  water 
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vapor  and  aerosols,  it  acts  as  an  excellent  absorber.   By 
equation  (13)  it  is  therefore  a  poor  reflector. 
3 .    Target  Reflectance 

The  characteristics  of  target  reflectance  follow  the 
same  behavior  as  the  wave  facet  reflectance.   The  difference 
in  the  magnitude  of  the  polarized  components  is  a  function 
of  the  index  of  refraction  of  the  surface  material  of  the 
target.   As  seen  in  previous  cases,  maximum  polarization 
effects  due  to  emission,  and  thus  polarization  will  occur 
nearer  to  grazing  incidence  angles.   The  net  polarization 
will  be  dependent  on  the  balance  between  thermal  emission 
and  reflection  of  incident  radiance,  related  through  Eq. 
(15) .   For  the  purpose  of  the  experiments  performed  to 
demonstrate  the  different  polarizing  imaging  techniques, 
target  distance  and  lack  of  detailed  geometric  orientation 
of  surfaces  precludes  exact  measurement  of  target  surface 
polarization;  however,   qualitative  effects  will  be  shown. 
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IV.   DATA  ACQUISITION  TECHNIQUES 

To  analyze  the  target  radiance  as  it  is  affected  by 
temperaturer,  polarization,  and  other  factors,  it  is 
necessary  to  employ  thermal  imaging  equipment  for  recording. 
The  split-field  polarimeter  and  internal  filters  use  the 
AGA-780  thermal  imaging  system  for  data  acquisition.   In 
the  case  of  the  internal  filters  which  were  designed  to  fit 
in  only  the  AGA-780  Thermovision,  the  choice  is  necessary, 
but  the  split-field  is  made  to  adapt  to  a  number  of 
different  sensors,  particularly  the  AMBER  4128  focal  plane 
array  system.   All  the  experimental  data  taken  to  be 
analyzed  was  taken  with  the  AGA-780  Thermovision. 
Accompanying  software  necessary  to  digitize  and  display 
radiance  data  consisted  of  the  AGA  CATS  software  and 
Interactive  Data  Language  (IDL)  programs.   The  following  is 
a  brief  discussion  of  the  data  collection  devices. 
A.    AGA-780  THERMOVISION 

The  AGA-708  Thermovision  consists  of  a  two  channel 
infrared  scanning  detector  system.   The  AGA  can  operate  in 

both  the  2-5.6(om  and  8-12(jm  bands  using  liquid  nitrogen 

cooled  InSb  and  HgCdTe  detectors,  respectively.   Scanning  is 
accomplished  by  two  rotating  prisms  which  scan  the  image 
across  the  detector,  yielding  1.1  mrad  geometric  resolution. 
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The  input  lenses  for  the  Long-wave  infrared  (LWIR)  used  for 
the  experiments  consist  of  anti-reflection  coated  Germanium. 
The  LWIR  band  was  used  for  all  imaging.   The  optical  system 
is  f/1.87,  with  a  7x7  degree  field  of  view.   This  f  number 
was  a  key  design  input  for  the  split  system  to  match  the 
aperture  of  the  AGA.   The  AGA-78  0  is  mounted  on  a  common 
plate  with  the  split-field  system  with  a  sliding  mount  to 
facilitate  easy  changing  between  the  LWIR  and  MWIR  band 
detectors.   Figure  11  in  chap.  7  shows  the  arrangement  used 
with  the  split  system. 
B.    AGA  CATS  SOFTWARE  PACKAGE 

The  analog  data  provided  by  the  scanning  detector  is 
digitized  for  display  and  measurement  by  the  AGA  CATS 
software.   The  CATS  software  uses  an  8  bit  A/D  converter  to 
convert   detector  intensity  into  a  16  color  display  for 
real-time  viewing  on  a  computer  screen.   Many  different 
functions  are  available  to  the  user  such  as  storing  frozen 
image  data,  image  manipulation,  temperature  readout  and 
emmissivity  compensation.   When  using  the  CATS  program  to 
provide  temperature  readings  of  a  target,  it  is  necessary  to 
calibrate  the  system  to  match  target  radiance  with 
temperature.   This  is  accomplished  by  recording  the  thermal 
level  seen  by  the  A/D  conversion  and  fitting  a  curve  to  the 
data  through  a  least-squares  method.   The  curve  is 


30 


reproduced  by  the   equation    [Ref   18] 


ThermalLevel— -  OS  (17 

C  eB/T-l 


A,    B,  and  C  are  determined  by  the  curve  fitting,   OS  is  the 
offset  correction  up  or  down,  and  T  is  the  temperature 
(ordinate)  of  the  plot.   Thermal  Level  value  range 
determines  the  abscissa.   Plots  for  system  calibration  for 
the  internal  filters  and  with  the  split-image  system  added 
are  shown  in  Appendix  D. 

The  most  important  feature  of  CATS  is  the  ability  to 
store  frozen  image  data  for  later  analysis.   To  make 
detailed  measurements  and  perform  data  manipulation,  further 
software  is  required  with  greater  flexibility. 
C.    INTERACTIVE  DATA  LANGUAGE  (IDL)  SOFTWARE 

In  order  to  process  the  image  data  to  a  level 
comensurate  with  the  accuracy  in  measurement  desired,  a 
software  package  is  required  that  provides  flexibility, 
digital  resolution,  and  data  manipulation  capability.   IDL 
allows  the  user  all  of  these  functions  and  is  required  to 
give  the  data  processing  power  needed.   It  is  produced  by 
Research  Systems  Incorporated,  Boulder,  CO.   Data  that  is 
stored  by  the  CATS  program  is  read  into  IDL  in  array  format 
where  numerical  analysis  is  performed  with  the  data 
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displayed  in  an  image  format.   Routines  are  written  to 
perform  a  number  of  functions  and  will  be  discussed  in  the 
analysis  chapter. 
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V.   POLARIZATION  FILTERING  TECHNIQUES 

Use  of  polarization  signature  and  background  features 
for  target  discrimination  implies  the  employment  of 
polarization  filters.   The  technology  for  large  aperture 
polarizers  in  the  IR  is  relatively  new.   Most  methods 
commonly  used  in  the  visible  range  are  not  available  in  the 
IR. 
A.    POLARIZING  FILTERS 

Polarizing  filters  have  traditionally  been  constructed 
using  anisotropic  crystals  or  wire  grids.   Of  these  methods, 
the  wire  grid  technique  has  been  applied  in  recent  years 
using  lithographic  techniques.   This  consists  of  a 
conducting  grid  deposited  on  a  transparent  substrate. 
Placed  in  the  optical  path  of  the  incoming  radiation,  the 
filter  will  pass  the  component  of  the  radiation  with  the 
electric  field  orthogonal  to  the  orientation  of  the  grid. 
The  grid  behaves  like  a  dielectric  sheet  in  this  case.   The 
component  of  the  radiation  with  the  electric  vector  parallel 
to  the  wire  grid  induces  currents  in  the  wires  and  is 
reflected.   Fig.  9  shows  a  schematic  of  this  process.   A 
restrictions  on  the  grid  is  that  the  line  spacing  must  be  of 
the  order  of  the  wavelength  of  the  incident  radiation 
(Ref.6).   Such  filters  for  the  infrared  are  available  with 
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clear  apertures  up  to  45mm.    Typical  transparent  substrates 
in  the  infra-red  are  materials  such  as  Germanium,  Silicon, 
and  KRS-5. 
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Figure  9.   Wire  Grid  Polarizer 

1.    External  Filters 

The  location  of  the  filter  in  the  optical  path  has  a 
large  impact  on  image  attributes.   If  a  flat  reflecting 
surface  is  placed  in  front  of  the  collimated  lens,  it  is 
possible  for  the  cold  detector  to  'see'  itself  in  the  image. 
This  is  referred  to  as  the  'Narcissus  Effect'  (Ref.23), 
which  perturbs  the  image  greatly.   Tilting  the  polarizer  may 
alleviate  this  problem  but  will  reflect  radiation  into  the 
image  from  outside  the  field  of  view.   Another  drawback  to 
using  an  external  filter  is  the  size  reguired  to  cover  the 
aperture  defined  by  the  input  lenses.   It  is  necessary  to 
cover  the  entire  aperture  which  is  a  function  of  the  system 
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f  number  and  focal  length.  Benefits  of  using  an  external 
filter  are  ease  of  insertion/removal  and  grid  orientation 
changing. 

The  large  external  filter  was  used  in  the  MAPTIP  exper- 
iment in  October  1993.   The  large  external  polarizer  pro- 
vided excellent  data  for  the  determination  of  the  degree  of 
polarization  of  the  targets  and  backgrounds  and  of  contrast 
improvement  (Ref.l).  Conclusions  showed  percent  polarization 
of  vertical  sea-water  emission  up  to  83%  in  the  long-wave  IR 
band.   Contrast  improvement  from  emission  polarization 
filtering  was  as  high  as  36.3%.   Large  improvements  such  as 
this  justify  the  utility  of  the  external  polarizers  and  show 
their  effectiveness. 

2.    Internal  Filters 

Filters  can  also  be  placed  inside  the  AGA-780  imaging 
system  between  the  main  input  lenses  and  the  detector.   This 
allows  hands-off  operation  and  smaller  filter  size.   In 
systems  which  are  configured  to  use  filters  in  the  optical 
path  for  other  types  of  filtering,  small  polarizing  filters 
can  be  substituted.   These  filters  will  typically  be  de- 
signed into  non-imaging  locations.   'Narcissus'  spots, 
filter  mounting  and  alignment  problems,  and  other  complica- 
tions can  thus  be  avoided  with  this  method.  Implementation 
can  be  accomplished  with  minor  system  modification  and  cost. 
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Installation  in  an  existing  filter  wheel  requires  a  polar- 
izer for  each  polarizer  orientation  to  be  used,  adding  to 
system  expense.     Internal  filtering  is  a  highly  desirable 
option  and  experimental  use  will  be  discussed  in  more  detail 
later. 

An  alternative  method  of  internal  polarizing  applicable 
to  focal  plane  array  systems  is  placing  an  independent 
filter  closely  against  each  element  of  the  focal  plane 
array [6] .   Each  element  of  the  array  can  then  have  any 
preselected  orientation  of  linear  polarization  built  into 
the  focal  plane  mask.   An  added  benefit  over  rotating 
polarizers  is  the  simultaneous  recording  of  all  the  selected 
polarizations  without  manipulating  a  filter  selector. 
Disadvantages  are  the  pixel  offset  of  the  image  for  each 
polarization,  and  the  inherent  loss  of  resolution.   The 

practicality  of  this  technique  for  the  3-5[im  band  has  been 
demonstrated  in  a  NASA- supported  program  (Ref .  6) . 
B.    SPLIT  FIELD  OPTICS 

A  method  of  polarization  in  which  the  polarizer  is 
internal  to  an  optical  system  but  external  to  the  imaging 
system  itself  is  the  technique  of  Split-Field  polarization 
filtering.   The  concept  of  the  split  field  system  is  that 
orthogonally  polarized  images  can  be  extracted  from  the 
field-of-view  simultaneously  with  one  polarizing  filter  and 
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presented  in  adjacent  sections  of  the  final  image  field  of 
view.   This  is  accomplished  by  splitting  the  beam  with  a 
polarizer.   The  Environmental  Research  Institute  of  Michigan 
(ERIM)  has  developed  such  a  system  which  uses  Brewster  angle 
reflection  to  separate  one  plane  polarized  component  from 
the  beam,  extracting  the  orthogonal  component  from  the 
transmitted  component  through  the  Brewster  plate  (Ref.  5). 
Since  the  Brewster  plate  is  tilted,  a  reference  is  needed  to 
keep  uncalibrated  background  from  contaminating  the  image, 
as  in  the  NPS  system.   Commonality  exists  between  the  ERIM 
system  and  the  NPS  system.  However,  the  NPS  system  defines 
the  apertures  of  the  two  channels  and  brings  them  to 
adjacent  halves  of  the  final  image  plane.  The  ERIM  system 
focusses  the  two  channels  onto  separated  areas  on  a  single 
integrated  focal  plane  array. 

The  Naval  Postgraduate  School,  Naval  Academic  Center 
for  Infrared  Technology  (NACIT)  has  constructed  and  begun 
testing  a  split-field  system.   Polarization  is  attained 
through  a  wire  grid  polarizer  tilted  to  the  incoming 
radiation  with  polarization  orientation  obstained  from  grid 
transmission  and  grid  reflection.   When  considering  options 
for  selecting  an  imaging  system  for  the  split-field  optics, 
the  scanning  AGA  was  chosen.   The  only  focal  plane  array 
available  for  use  at  this  time  operates  in  the  3-5  pirn  band. 
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Because  of  the  nature  of  the  targets  and  environments  where 
the  split  system  will  be  used  most,  it  was  decided  to  first 
build  the  system  to  operate  in  the  8-12  \xm   region.   The 
existing  hardware  and  familiarity  using  the  AGA  coupled  with 
IDL  software  and  established  data  analysis  techniques  made 
it  a  good  choice  for  examining  the  new  system.   The  trade  to 
go  to  the  AGA  is  in  reduced  frame  sampling  rate  and 
resolution.   A  detailed  explanation  of  the  split-field 
technique  employed,  a  system  description,  and  polarizing 
characteristics  of  the  system  will  follow  later. 
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VI 


INTERNAL  POLARIZING  FILTERING  WITH  THE  AGA-780 


Internal  filtering  of  the  radiation  after  it  passes 
through  the  imaging  system  input  lenses  can  be  accomplished 
with  the  AGA-780  Thermovision  Scanning  Imaging  system  by 
using  an  existing  filter  wheel.   The  original  purpose  of  the 
filter  wheel  was  to  hold  various  accessory  spectral  filters 
in  the  optical  path  before  the  radiation  reached  the 
detector.   Fig.  10  shows  the  internal  layout  of  the  AGA-780. 
One  of  seven  positions  of  the  wheel  is  selected  by  rotating 
the  knob  on  the  front  face  of  the  unit.   Care  must  be  taken 
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Figure  10.   Internal  Layout  of  the  AGA 
Thermovision  (model  with  only  one  IR  Band) 
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to  ensure  that  no  internal  circuitry  compensates  for  an 
assumed  filter  present  at  the  position  of  the  polarizer. 
Polarization  measurements  with  this  system  require 
sequential  images  with  alternate  filter  wheel  positions 
inserting  vertically  and  horizontally  polarized  filters. 
This  necessarily  involves  time  delays. 
A.    ALUMINUM/KRS-5  FILTERS 

The  filters  used  consisted  of  an  aluminum  grid 
superimposed  on  a  9.5mm  diameter  x  3mm  thick  round  KRS-5 
substrate  manufactured  by  Graseby-Specac;  Suffolk,  England. 
The  performance  specifications  provided  with  the  filters  are 
nearly  identical  for  both  and  will  be  treated  as  such  for 
explanation.   Appendix  A  shows  the  performance  curves  for 
the  two  filters.   The  various  traces  numbered  1-3  represent 
power  transmitted  with  the  grid  vertical,  horizontal,  and 
with  the  grids  from  two  filters  crossed,  relative  to  the 
electric  vector  of  the  incoming  radiation.   The  determina- 
tion of  the  perpendicular  and  parallel  transmittances  for 
the  filter  involves  use  of  a  spectrometer  which  itself 
includes  internal  polarizing  reflections.   Thus  it  must  be 
accomplished  by  measuring  the  transmittance  twice,  with  the 
grid  in  two  orthogonal  orientations.   These  two  measured 
transmittances  are  designated  Ev  and  Eh.   The  corresponding 
instrument  transmittance  is  measured  without  the  grid  in 
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place  as  E0,  and  the  transmittance  for  two  crossed  grids  as 
E1.       From  these  quantities  the  grid  transmittances  are 
derived  as 
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These  quantities  are  specified  by  the  manufacturer  for  each 
filter.   The  expressions  are  approximate  values  derived  by 
Young,  et  al .  (Ref .  13)  . 
B.    FILTER  INSTALLATION  IN  THE  AGA-780 

Figure  10  shows  the  internal  layout  of  the  AGA-780.   To 
install  the  filters,  the  appropriate  holes  in  the  wheel  were 
enlarged  approximately  0.5mm  in  diameter.   Two  filters  were 
then  installed  in  filter  positions  4  and  5,  carefully 
oriented  with  passage  axis  vertical  and  horizontal, 
respectively.  The  filter  wheel  was  then  installed  in 
accordance  with  the  manual.   Extreme  care  was  taken  to 
select  positions  in  the  filter  wheel  which  have  no 
electronic  filter  compensation  in  the  system.   The  manual 
denotes  which  positions  these  are.   Once  the  system 
modification  was  complete,  new  calibration  constants  were 
computed  and  installed  in  the  CATS  system.   Calibration 
constants  are  listed  in  Appendix  D. 
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VII.   THE  NPS-NACIT  SPLIT-FIELD  POLARIMETER 

The  split-field  system  approach  possesses  many 
advantages  over  external  and  internal  filter  placement 
normal  to  the  path  of  the  radiation,  particularly  in 
elimination  of  time  lags  and  registration  problems  between 
orthogonal  polarizations.   The  Naval  Academic  Center  for 
Infrared  Technology  (NACIT)  at  the  Naval  Postgraduate  School 
(NPS)  has  constructed  and  tested  a  split-field  front  end 
optical  system  for  input  to  thermal  imaging  systems, 
specifically  matched  to  the  AGA-780  and  the  AMBER 
Engineering-4128 .   The  Split-Field  Polarimeter  is  the  result 
of  a  conceptual  design  by  Professor  Alfred  W.  Cooper.   The 
actual  construction  and  design  optimization  was  performed  by 
Professor  Eugene  C.  Crittenden  and  LT  Peter  H.  Heisey,  USN, 
as  assistant.   The  general  features  of  the  NPS-NACIT  Split- 
Field  Polarimeter  are: 

1.  Simultaneous  display  and  recording  of  orthogonally 
polarized  images. 

2.  Utilization  of  entire  f ield-of-view  for  display. 

3.  Adaptability  to  different  imaging  systems. 

4.  Dual-band  (8-12(xm  and  3-5|im)  operation  using 
interchangeable  filters. 

5.  F/1.87,  99mm  focal  length  optics. 


43 


As  mentioned  in  Chapt .  IV  for  other  systems  with  split-field 
optics,  the  'Narcissus  spot'  is  avoided  by  having  the  filter 
tilted  relative  to  the  path  of  the  radiation.   A  background 
plate  is  also  needed  to  give  a  uniform  reflected  radiation 
source  into  the  images.   The  specifics  of  the  system  and  its 
components  and  design  will  be  discussed. 
A.    SYSTEM  DESIGN 

The  split-field  optical  system  consists  of  input  and 
output  compound  lenses,  polarizing  splitter,  and  front- 
surface  aperture  limiting  mirrors.  The  components  are 
arranged  in  a  vertical  plane  with  respect  to  the  ground. 
The  system  takes  the  collimated  input  radiation,  focuses  it 
down  as  it  is  split  into  two  orthogonally  polarized 
components,  and  folds  and  recombines  the  two  split  halves  of 
the  image  while  recollimating  the  light  for  processing  by 
the  thermal  imaging  system.   In  this  process,  the  optical 
aperture  is  redefined  so  that  only  the  scan  lines  of  the 
central  half  of  of  each  image  are  retained.  The  remaining 
3.5x7  degree  halves  are  combined  in  adjacent  halves  of  an 
intermediate  7x7  degree  image  plane.   The  7x7  degree 
field  of  view  of  the  optical  system  is  thus  split  into  two 
channels  with  orthogonal  polarizations.   Each  channel  is 
aperture  limited  into  a  7x3.5  degree  field  to  compose  half 
of  the  final  image.   The  original  radiant  flux  in  the  image 
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is  preserved  in  this  way.   An  overall  picture  of  the  system 
in  use  with  the  AGA-780  Thermal  Imaging  System  is  shown  in 
Fig.  11.   In  this  setup,  the  entire  imaging  system  is 
mounted  on  one  lightweight  mounting  plate.   Fig.  12  shows 
the  inside  of  the  box  with  all  components  exposed.   Aluminum 
was  used  for  all  mountings  and  plates  because  of  low  weight, 
corrosion  resistance,  and  easy  machining. 
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Figure  11.  Top  View  of  AGA  with  Split-System. 
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Figure  12.   Side  View  of  AGA  with  Split-System, 


Figure  13.   Front  View  of  Split-Field  input 
lenses  with  AGA  in  background. 
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availability,  corrosion  resistance,  and  easy  machining.   The 
front  view  in  Fig.  13  shows  the  exposed  input  lens.   The 
long  wave  input  lens  for  the  AGA  system  is  hidden  by  the 
backside  of  the  system.   The  short  wave  lens  not  being  used 
by  the  AGA  is  to  the  side  and  may  be  slid  into  place  behind 
the  optics  when  a  short-wave  polarizer  is  installed.   Figure 
14  shows  the  image  orientation  and  focal  plane  within  the 
split-field  optics. 
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Figure  14.  Schematic  of  Ray  Paths  in  Split  Optical  System 

The  incoming  rays  encounter  two  Germanium  lenses  (labeled 
lens  system)  and  are  refracted  to  a  focal  plane.   In  the 
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process,  the  polarizing  splitter  splits  the  beam  into  two 
paths  with  one  beam  vertically  polarized  with  respect  to  the 
ground  (transmitted  beam)  and  the  other  reflected  with 
horizontal  polarization.   The  transmitted  beam  is  reflected 
to  a  focal  plane  at  the  edge  of  a  small  mirror  used  by  the 
reflected  beam  to  be  turned  in  the  same  direction  (Fig.  14) . 
At  this  point,  both  beams  have  been  recombined  in  a  single 
focal  plane.  The  beams  then  diverge  and  are  reflected  by  the 
last  mirror  into  two  lenses  identical  to  the  first  lens 
system  (labeled  lens  system) .    The  radiation  emerges 
afocally  so  that  the  AGA  imager,  positioned  immediately 
behind  the  system  can  process  the  image  for  display. 
B.    SPLIT-FIELD  OPTICAL  COMPONENTS 

The  split-field  system  consists  of  3  major  components 
as  stated  in  the  previous  chapter.   The  specific 
characteristics  and  operation  of  these  components  will  now 
be  addressed. 

1 .    Lenses 

The  lenses  used  by  the  system  to  generate  the  split 
image  for  input  to  the  thermal  imaging  system  consist  of  two 
identical  assemblies  each  containing  two  germanium  anti- 
reflection  coated  lenses.  One  assembly  is  used  for  the  input 
and  the  other  for  the  output  side.  The  output  assembly  is 
reversed,  end  for  end,   relative  to  the  input  assembly. 
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The  inner  surface  of  the  first  lens  of  the  input  assembly 
and  the  inner  surface  of  the  last  lens  of  the  output 
assembly  are  aspheric.  This  reduces  aberrations  and  provides 
a  minimum  spot  size  of  1.1  mrad,  within  the  resolution  of 
the  AGA  imager.  The  lens  surfaces  are  coated  for 

antiref lection  in  both  the  4|im  and  12|im  wave  bands.   This  is 
necessary  because  the  high  index  of  refraction  of  germanium 
(4.0)  would  produce  large  reflection  losses  and  scattered 
radiation  if  the  surfaces  were  left  uncoated.  The  overall 
focal  length  of  each  of  the  two  lens  assemblies  is  99.2mm 
with  an  f  number  of  1.87 

2.    Polarizing  Splitter 

The  most  important  component  of  the  split  system  from 
the  polarizing  standpoint  is  the  Al/Ge  polarizing  splitter. 
The  custom  designed  splitter  consists  of  a  germanium 
substrate  with  an  aluminum  grid  deposited  on  one  side.   The 
grid  consists  of  .25(im  wide  aluminum  strips  separated  by 

.40nm  for  an  overall  period  of  .  65|im.   The  substrate  is  a 
round  disc,  50mm  in  diameter  x  2mm  thick.  Both  sides  of  the 
splitter  areanti-ref lection  coated  for  the  8-12|im  band.   It 

is  necessary  to  replace  the  splitter  for  the  short-wave  IR 
band,  since  the  vendor  found  it  impossible  to  coat  the 
element  for  both  bands.   The  same  analysis  done  for  the  KRS- 
5  polarizers  in  Chapt.  V  was  performed  for  the  Ge  splitter. 
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Data  provided  by  the  manufacturer  (Graseby-Specac)  is  valid 
for  normal  incidence  of  the  optical  beam  on  the  splitter. 
Values  for  Klr    K2,  and  P  are  .81,  0.0,  and  1.0,  respectively 
at  10|im.   This  means  the  splitter  passes  81%  of  the  linearly 
polarized  light  perpendicular  to  the  grid  orientation  and  is 
~100%  polarized.   The  manufacturers  performance  curves  are 
provided  in  Appendix  A.   Extensive  analysis  was  done  to 
verify  operation  of  the  splitter;  this  will  be  addressed  in 
section  C.   Results  showed  that  there  is  a  wavelength 
dependency  of  the  radiation  fraction  that  is  reflected  and 

transmitted.   It  is  found  that  wavelengths  between  7|im  and 
9|im  favor  reflection  and  the  opposite  is  true  above  9fjm. 

This  is  due  to  inefficiencies  in  the  polarizing  grid  and 
anti-reflection  coating.    Results  obtained  by  system 
calibration  show  this  and  will  be  discussed  later. 

3.    Gold  Mirrors 

The  folding  operation  of  the  optical  system  is 
performed  by  four  gold  first-surface  'protected'  mirrors. 
The  mirrors  were  purchased  from  Edmund  Scientific  Co.  and 
are  advertised  to  be  1/10  X   flat  at  .5 pun.   Gold  was  chosen 
because  of  its  high  reflectivity,  and  thus,  low  emissivity. 
Figure  15  shows  a  plot  of  the  emissivities  at  9.92  urn  for 
gold  at  a  range  of  incidence  angles. 
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Figure  15.   Emissivity  of  Gold  from  eq. (9) . 

As  expected,  significant  polarization  effects  occur  only  at 
angles  close  to  grazing  incidence.   Most  reflecting  angles 
of  the  mirrors  are  near  45  degrees,  so  these  effects  are 
small  (<1%) ,    but  must  be  considered. 
C.    SYSTEM  CHARACTERISTICS 

The  combined  effects  of  the  presence  of  the  components 
in  the  optical  path  and  mechanisms  by  which  they  perform 
their  functions  must  be  quantified  to  predict  overall  system 
performance.   Irregularities  in  refracting  and  reflecting 
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surfaces,  polarized  reflection  and  emission,  and  system 
background  radiance  are  effects  which  influence  the  measured 
radiance  and  resolution  of  the  target.   Experiments  have 
been  performed  to  characterize  these  effects  to  assist  in 
predicting  final  image  quality. 

1 .    Resolution 

To  obtain  the  full  resolution  possible  by  the  AGA-780 
when  used  with  the  split-field  system,  effects  induced  by 
imperfections  in  the  optics  of  the  split-system  must  be 
minimized.   Upon  delivery  of  the  splitter,  visual  resolution 
was  tested  using  a  standard  bar  chart.   Reflection  of  the 
bar  chart  off  the  surfaces  of  the  splitter  showed  no  loss  in 
resolution  from  direct  viewing  of  the  chart.   It  follows 
that  imperfections  in  the  surfaces  would  induce  less  effect 
in  the  longer  wavelengths  and  would  not  be  considered. 

A  test  was  also  done  on  the  combined  AGA-split  system 
to  discern  if  resolution  was  degraded  by  the  presence  of  the 
split-field  system.   A  sharp  edged  (cold)  plate  was  placed 
in  front  of  a  black-body  thermal  source.   Using  the  AGA 
optics  alone,  measurements  were  taken  of  the  number  of 
pixels  in  the  image  that  were  required  to  resolve  the  sharp 
edge  between  hot  and  cold  surfaces.   Along  the  direction  of 
image  raster  scan,  the  resolution  is  about  2-3  pixels. 
Perpemdicular  to  the  raster  scan  the  resolution  is  slightly 
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better.  Parallel  to  the  raster  scan  the  resolution  is 
apparently  limited  by  the  frequency  band  pass  of  the 
amplifier  electronics  and  the  response  time  of  the  detector. 
In  the  direction  perpendicular  to  the  raster  scan,  the 
resolution  is  slightly  higher  because  the  response  times  of 
the  detector  and  amplifier  are  not  involved.  In  both 
parallel  and  perpendicular  scan  directions,  no  loss  of 
resolution  occurred  in  the  image  when  the  split-field  optics 
were  added. 

2.  MTF 

Modulation  transfer  function  tests  were  performed  for 
the  AGA  system  with  and  without  the  split-field.  The  cutoff 
frequency  for  the  AGA  in  the  direction  of  raster  scan  is  . 3- 
.5  cycles/mrad  which  gives  a  resolution  of  2-3mrad.   The  AGA 
system  manual  gives  a  system  resolution  of  l.lmrad.  This  is 
apparently  the  geometrical  resolution  of  the  AGA  optical 
system.  The  MTF  data  show  that  the  overall  system  has  a  2-3 
mrad  resolution.  There  was  no  loss  in  spatial  resolution 
with  the  split-field  system; 

3.  System  Polarization  Performance 

The  main  driver  in  the  performance  of  the  system  is  the 
polarizing  splitter.   Establishing  numerical  values  for 
transmissivity,  reflectivity,  and  absorption  of  the  splitter 
for  polarized  and  unpolarized  target  radiation  requires 
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measurement  of  all  radiance  sources  throughout  the  image  and 
split-field  system. 

a.         Spectral    Polarizer  Performance 
The  manufacturer  of  the  polarizer  provided 
spectral  transmittance  data  for  the  polarizer  at  normal 
incidence  to  a  reference  beam.   For  use  in  the  split  system, 
it  is  necessary  to  know  the  spectral  transmittance  and 
reflection  for  45  degree  incidence. 

To  measure  the  transmittance,  the  polarizer  was  placed 
at  45  degrees  to  the  incoming  beam  in  a  Perkin-Elmer  Model 
337  Infrared  Grating  Spectrophotometer.   The  photometer 
measures  the  characteristics  in  the  same  manner  as  the 
spectrophotometer  used  by  the  manufacturer.   Because  of 
polarization  induced  in  the  beam  by  the  grating,  the  method 
described  by  Young  et  al . (Ref  13) ,  is  used  to  compute  the 
transmittance.   This  is  the  approach  outlined  in  Chapt .  VI. 
Appendix  B  shows  the  plots  generated  by  the  polarizer  at 
normal  and  45  degree  incidence.   Agreement  with  the 
manufacturers  values  are  within  3%  for  normal  incidence. 
For  wavelengths  less  than  10.5|am,  the  transmissivity  is 
higher  than  at  normal  incidence.   This  is  probably  due  to 
the  fact  that  the  effective  grid  spacing  is  reduced  by 
cos  (45)  ,  when  the  grid  is  perpendicular  to  the  plane  of 
incidence  of  the  beam,  thus  moving  the  curve  left  to  shorter 
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wavelengths.   It  is  necessary  to  mention  that  placing  the 
filter  in  such  a  manner  as  to  reflect  ambient  radiation  into 
the  detector  does  not  contaminate  the  readings.    The 
intensity  radiated  by  a  typical  room  temperature  object 
compared  to  that  of  the  IR  source  in  the  device  is  on  the 
order  of  0.001. 

Obtaining  performance  for  reflection  polarization 
required  altering  the  conventional  configuration  of  the 
polarizer  to  the  beam  to  utilize  the  front  surface  as  the 
effective  transmitter.   Figure  16  shows  the  arrangement  used 
to  obtain  this. 
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Figure  16.   Configuration  for  polarization 
reflection  measurement  using  the  Perkin-Elmer 
M337  Infrared  Grating  Spectrophotometer. 
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Data  for  the  linear  polarization  reflection  efficiency  is 
provided  in  Appendix  B.   Reference  13  provides  a  derivation 
of  the  measurement  technique  to  obtain  the  transmittance  of 
a  polarizing  grid.   The  difference  in  finding  the  efficiency 
here  is  that  reflection  is  substituted  for  transmission  in 
the  equations.   The  radiant  intensity  with  the  polarizer 
removed  was  obtained  by  placing  an  aluminum  mirror  in  place 
of  the  polarizer  and  resetting  the  100%  transmission 
indicator.   The  performance  plot  in  appendix  B  shows  this 
trace  across  the  top  of  the  graph.   The  designation  for  the 
reflected  efficiency  is  K1R.   The  equations  are  analogous  to 
those  described  in  chapter  VI.   It  was  not  possible  to 
measure  K2R  (orthogonal  polarization  reflection)  because  of 
the  unique  arrangement  and  nonavailability  of  another 
similar  polarizer. 

The  degree  of  reflected  polarization  perpendicular 
to  the  grid,  can  be  described  as  follows.   Any  radiation 
impinging  on  the  grid  with  E  vector  perpendicular  that  is 
not  transmitted  must  be  reflected  due  to  fact  that  the  grid 
width  is  small  compared  to  the  wavelength,  precluding 
current  generation.   This  reflected  component  is  the 
complement  to  the  transmitted  fraction  of  the  radiation. 
This  theory  will  be  demonstrated  in  the  following  section. 
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b.         Polarizer  Performance   Utilizing  Radiance 
Analysis 

(1)  Theoretical  Representation 

An  alternative  method  for  measuring  the  fraction 
of  polarized  radiation  transmitted  and  reflected  by  the 
polarizing  splitter  is  to  take  radiance  readings  of  the 
radiation  passed  with  a  known  source  polarization 
orientation.   However,  complications  arise  with  this  method 
when  target  temperature  is  near  the  temperature  of  the 
optical  system,  because  the  background  irradiance  that 
enters  the  detector  is  of  the  order  of  the  target 
irradiance.   To  discriminate  the  target  irradiance  from  the 
background,  a  mathematical  description  of  the  background 
effect  is  required. 

For  the  case  of  unpolarized  radiation  incident  on 
the  optical  system,  the  following  equations  apply: 


J        =f„  I+f-  J+S„ 

*,      I  R±    x         Rt 


(19) 


R 


VNPOL  "1       *  ■!■  B 


where 


S  -  ,  £_ —the  respective  background  emittances 

TB  RB 

measured  for  each  channel,    independent 
of  target  irradiance 
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e    ,  s    =e        +g  +g  +g 

TD         RD  lenses       mirrors       splitter       background  plate 


lB         ^B 


X^and  X.  represent  the  components  of  target  irradiance 

perpendicular  and  parallel  to  plane  of  incidence  of  the 
polarizing  splitter.   These  are  what  is  actually  measured  by 
the  thermal  imaging  system.   For  the  entire  analysis  the 
convention  used  is  that  the  polarizing  splitter  grid  lines 

are  perpendicular  to  the  plane  of  incidence.   The  eTB,RB  terms 

represent  the  radiant  exitance  of  all  the  components  of  the 
system  when  the  incident  intensity  into  the  system  ~0.   This 
is  simulated  by  directing  the  optics  at  a  pool  of  liquid 
nitrogen.   The  subscripts  T  and  R  designate  the  transmitted 
and  reflected  channel  through  the  split-field  optical 
system.   T  denotes  the  radiance  fraction  which  passes 
through  the  polarizing  splitter  and  is  displayed  as  the  top 
half  of  the  split  image;  R  denotes  the  other  channel 
reflected  by  the  splitter.   The  perpendicular  and  parallel 
subscripts  describe  the  orientation  of  the  E  vector  for  the 
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irradiance  relative  to  the  plane  of  incidence  of  the 
splitter.   For  example,  unpolarized  radiation  is  commonly 
represented  by 


t^r^^J^+^^J2) 


20 


To  eliminate  unknown  terms  in  Eq.  (19) ,  measuring 
a  purely  linearly  polarized  source  in  two  orthogonal 
orientations  yields  a  simpler  set  of  equations.   For  the 
first  case  of  purely  polarized  target  irradiance 
perpendicular  to  the  plane  of  incidence  (perpendicular  to 
the  grid) , 


X  =£„  J, +6+8, 

r«   r'  '   r*   F'i  (2i 


The  fT/R  terms  denote  that  fraction  which  goes  into  the 
respective  channels  from  splitter  reflection  and 
transmission  performance.   The  target  consisted  of  a 
polarizing  filter  in  front  of  a  black  body  source; 
therefore,  there  is  some  filter  emittance  that  must  be 

considered.   The  added  ef  terms  represent  the  emittance  of 
the  polarizing  source  for  each  channel,  for  each  orientation 


59 


of  the  target  producing  filter.  It  can  be  seen  that  there  is 
a  fractional  contribution  to  each  channel  to  be  measured   T 
and  R.   This  is  because  it  is  known  that  the  transmittance 
of  the  grid  is  not  1.0.   For  the  case  where  the  target 
radiance  is  plane  polarized  perpendicular  to  the  plane  of 
incidence  (parallel  to  the  grid) , 


**  ^    '*  (22 

i=f„  r+e.  +s 


There  are  few  terms  seen  in  the  T  channel  because  it  is 
known  from  the  spectral  performance  curves  that  the 
transmitted  E  components  are  ~100%  polarized.  Equation  21 
shows  the  intensities  seen  in  the  two  channels  with  a 
vertically,  parallel  (p)  linearly  polarized  source  and  eq. 
22  with  90  degree  rotation  of  the  target  E  vector,  s 
polarized. 

(2)   Experimental  Radiance  Measurement 
To  compute  the  performance  of  the  system  using  the 
equations  in  the  previous  section,  it  is  necessary  to  take 
careful  radiance  measurements  of  each  term  in  the  equations 
with  a  linearly  polarized  target.   Appendix  C  contains  the 
values  obtained  in  the  laboratory  for  these  measurements. 
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To  compare  the  performance  of  the  polarizing  splitter 
channels,  the  values  in  Appendix  C  are  plugged  into  eqns 
(21)  and  (22)  and  then  subtract  the  equations  from  each 
other.   The  relationships  obtained  are 


<VVTi=21-9 


f   1=22.4 


Measurements  showed  that  there  is  no  dependency  of  target 
intensity  on  polarization  orientation,  so  I.~£    • 
Therefore, 

T'  *'  -  =  0.98  (23 


f*. 


This  equation  allows  calculation  of  the  reflected  fraction 
of  the  incident  radiation  which  is  normally  passed  through 
by  the  grid.   The  implications  of  Eq.  (23)  are  that  the 
transmittance  of  linearly  polarized  irradiance  perpendicular 
to  the  grid  is  always  greater  than  the  reflection  efficiency 
of  the  orthogonal  component  from  the  grid  front.   This  is 
seen  with  the  performance  curves  in  appendix  B.   It  also 
allows  correlation  of  splitter  performance  using  the 
radiance  method  with  that  seen  in  the  spectral  performance 
curves.   To  compare  the  radiance  method,  it  will  be 
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assumed  that  the  non-transmitted  E  vector  perpendicular  to 
the  grid  undergoes  Snell  reflection  and  fT+fR=l.   The 
spectrophotometer  data  measured  in  the  previous  section 

shows  that  at  9.52um,  the  above  expression  =  0.97. 
Agreement  is  seen  between  the  two  methods  at  that 
wavelength,  but  the  ratio  varies  strongly  at  other 
wavelengths.   Schematic  representation  of  what  happens  at 
the  splitter  within  the  optical  system  is  shown  in  Figs.  17- 
19.   Figure  17  indicates  the  fractional  performance  in  the 
LWIR  for  unpolarized  input.   For  this  case,  58%  of  the 
background  is  reflected  into  the  transmitted  (T)  channel  and 
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Figure    17.       Schematic   of   polarizing   splitter 
for  unpolarized   incident   radiation.    Grid 
perpendicular    to   plane   of   incidence. 
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combined  with  41%  of  the  target  irradiance.  The  reflected 
(R)  channel  combines  42%  of  the  background  with  59%  of  the 
target  irradiance.   Figure  18  indicates  the  fractional 
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Figure  18.   Schematic  of  polarizing  splitter 
for  horizontally  polarized  incident 
radiation.  Grid  perpendicular  to  plane  of 
incidence . 


performance  in  the  LWIR  for  horizontally  polarized  input. 
For  this  case  the  transmitted  beam  contains  58%  of  the 
irradiance  from  the  background  plate  and  essentially  no 
target  radiation.   Figure  19  shows  the  equivalent  situation 
for  vertically  polarized  input,  which  shows  83%  of  the 
target  irradiance  with  58%  of  the  background  irradiance  in 
the  transmitted  (T)  beam,  and  17%  of  the  target  irradiance 
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Figure  19.  Schematic  of  polarizing  splitter 
for  vertically  polarized  incident  radiation. 
Grid  perpendicular  to  plane  of  incidence. 


with  42%  of  the  background  plate  irradiance  in  the  reflected 
(R)  beam. 

The  radiance  analysis  of  the  polarizing  splitter 
performance  provides  a  gross  analysis  and  verification  of 
the  data  provided  by  the  manufacturer  and  by  the  grating 
spectrophotometer.   This  representation  of  transmittted  and 
reflected  target  irradiance  must  be  considered  when 
processing  split-field  images. 
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VIII.   IMAGE  ANALYSIS 

To  prove  the  concepts  presented  and  to  demonstrate  the 
effects  of  polarized  filtering  on  image  contrast,  it  is 
necessary  to  perform  a  preliminary  image  analysis  on  images 
recorded  with  the  thermal  imaging  system  in  two 
configurations;  internal  polarizing  and  split-image  optical 
system  polarizing.  Two  separate  experiments  were  conducted 
to  demonstrate  system  performance. 

In  March  1995,  1338  images  were  taken  during  a  field 

experiment  with  the  R/V  POINT  SUR  from  the  Monterey  Bay 

Aquarium  Research  Institute  marine  operations  building  at 

Moss  Landing,  CA.   All  of  the  polarization  filtering  was 

done  internally.   In  August  of  the  same  year,  images  were 

recorded  using  the  NACIT  Split-field  Polarimeter  from  the 

Naval  Postgraduate  School  beach  property  of  the  USS  John 

McCain,  an  Arleigh  Burke  Class  DDG.   The  images  recorded 

show  contrast  improvement  effects  from  polarization 

filtering. 

A.    DEFINITION  OF  CONTRAST  AND  TARGET -BACKGROUND  RADIANCE 
DIFFERENCE 

When  an  observer  views  an  image,  the  critical  parameter 

for  detecting  and  recognizing  the  target  amidst  the 

background  is  target  contrast.   Numerous  definitions  are 
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currently  used  to   define   contrast.      Common   examples   of   these 
are    [23] 

N         -N  N         -N 

Contrast=     target    gffggggg;     Contrast^     tarqet    fcac^roand 

H  (24) 

background  fJ  "t"7\f 

target        background 

Here  N  represents  radiance.   It  is  desirable  to  compare  the 
contrast  of  images  which  may  be  separated  in  time.   In  this 
case,  changes  in  system  gain  cause  the  apparent  digitized 
image  values  to  change  for  similar  viewing  conditions. 
Also,  changing  instrument  settings  changes  N  for  every  image 
when  the  controls  are  altered.   Therefore,  when  target  and 
background  conditions  vary  little  between  images,  the 
difference  in  N  may  be  a  more  accurate  indicator  of 
contrast.   Differences  in  filter  transmittance  cause  the 
conventional  contrast  equations  to  be  in  error  so  that  any 
artificial  thermal  level  changes  induced  by  this  are  ignored 
for  filters  with  the  same  calibration  curve  slope.   This  is 
the  case  for  the  internal  filters  (Appendix  E) .  To 
accommodate  this,  the  analysis  will  only  consider 
differences  in  the  apparent  radiances  of  the  target  and 
background.   The  contrast  parameter  to  be  measured  is  termed 
Target-Background  Radiance  Difference  (TBRD) .   This  is  the 
numerator  of  both  contrast  equations  in  Eq.  (24) . 
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TBRD  =  AN  -  <N+         >   ""    <Nh     .  >  (25; 

target  background 


The  average  values  of  N  are  used  because  of  the  varying 
composition  of  the  image  with  target  distance,  system  noise, 
and  measurement  points  in  the  image.   If  the  image  radiation 
is  polarization  filtered,  the  resultant  TBRD  is  a  function 
of  the  amount  of  filtering  only  if  all  optical  components 
and  scene  characteristics  remain  constant.   It  has  been 
shown  that  many  different  components  of  an  image  can  display 
different  degrees  of  polarization.   If  background  irradiance 
is  suitably  filtered  more  than  target  irradiance,  TBRD  will 
increase. 

When  comparing  two  images,  it  is  useful  to  quantify  the 
increase  in  image  contrast  improvement.   The  TBRD 
improvement  factor  is  defined  by 

AK 
TBRD.  =AN.  = i  (2  6) 

improvemen  t  imp       A  AT 

2 

B.    FIELD  MEASUREMENTS  WITH  INTERNAL  FILTERS 

From  May  16-23,  1995,  polarized  thermal  imaging 
experiments  were  conducted  from  the  Monterey  Bay  Aquarium 
Research  Institute  using  the  R/V  POINT  SUR  as  target. 
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1338  images  were  recorded  using  internally  filtered  and 
unfiltered  images  in  the  8-12um  wave  band.   The  ship 
presented  all  aspects  with  ranges  varying  from  .25  nmi  to 
the  range  limit  of  the  day.   Appendix  E  carries  a  complete 
listing  of  the  images  and  a  small  picture  of  each  with  the 
filename  shown  on  the  bottom  of  each  image.   Vertically  and 
horizontally  polarized  images  will  be  compared  for  target- 
background  radiance  difference. 

1.    Visual  Image  Comparison 

Selected  pairs  of  images  demonstrate  the  visual  TBRD 
and  contrast  improvement  between  vertically  and  horizontally 
polarized  images.   Unpolarized  images  were  not  considered 
for  this  comparison  because  of  the  nature  of  the  data 
analysis  process,  and  the  change  in  system  settings 
necessary  to  obtain  useable  images  due  to  the  changing  of 
system  gain  and  altering  the  byte  level  of  the  images 
significantly.   The  definition  of  TBRD  would  not  be  valid 
for  comparison  to  the  polarized  images  and  further 
compensation  to  TBRD  would  be  necessary.   The  images  are 
selected  from  the  compilation  listing  in  Appendix  E.   The 
contrast  between  the  two  orthogonally  polarized  members  of 
each  pair  of  images  in  Fig.  20a  and  20b  shows  the  lower 
apparent  radiance  of  the  sea  water  background  in  the 
horizontal  polarization.   It  is  necessary  to  mention  that 
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the  images  were  printed  on  a  Hewlett-Packard  LaserJet  4M 
600dpi  printer.   Differences  in  the  images  are  more  easily 
seen  in  a  grayscale  color  table  on  an  SVGA  monitor.   The 
bright  area  seen  in  the  top  right  corner  of  each  image  is  a 
constant  temperature  source  created  by  inserting  a  mirror  in 
the  field  of  view,  reflecting  a  calibrated  black  body 
source.   This  provides  an  unpolarized  reference  in  each 
image  to  compare  filter  transmittance  and  provide 
temperature  calibration  confirmation.   The  source  is  not 
discussed  in  this  analysis. 
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Figure  20a.   Comparison  of  vertically  and  horizontally 
polarized  LWIR  Band  (8-12jnm)  images  of  RV  POINT  SUR. 


70 


M5L23014  23May95 

Polarization:  Horizontal 

TBRD=22.8 


M5L23015  23May95 

Polarization:  Vertical 
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Polarization:  Vertical 
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Figure  20b.  Comparison  of  vertically  and  horizontally 
polarized  LWIR  Band  (8-12um)  images  of  R/V  POINT  SUR 
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2.    Numerical  TBRD  Measurement. 

The  images  presented  in  Fig.  20  were  extracted  groups 
of  images  with  the  same  aspects,  ranges,  and  path 
conditions.   The  average  radiances  are  calculated  by 
averaging  over  the  target  and  background  areas  within  an 
area  around  the  ship.   Appendix  F  contains  the  IDL  code  to 
perform  this  operation.   To  compare  the  TBRD  for  a  number  of 
images,  the  images  must  be  taken  with  target  and  background 
conditions  as  constant  as  possible  for  trends  in  the  data  to 
appear.   Figures  21-23  show  TBRD  measurements  made  on 
different  sets  of  images  grouped  by  ship  range  and  aspect. 
The  vertical  axis  gives  the  TBRD  in  dimensionless  units. 
The  8  bit  A/D  converter  sets  the  maximum  thermal  level 
reading  based  on  system  gain  to  255.   This  irradiance 
reading  is  independent  of  system  calibration  settings  and 
measures  the  raw  data  from  the  detector  with  no  compensation 
for  temperature  accuracy.   The  horizontal  axis  gives  range 
computed  from  ship  image  height  measured  in  pixels.  Smaller 
targets  appear  to  the  right  side  of  the  plot,  which  shows 
TBRD  as  a  function  of  range.   The  general  trend  of  TBRD 
reduction  with  target  size  is  due  to  increased  path 
attenuation  and  path  radiance  with  range.  The  legend  on  the 
plots  denotes  the  polarization  and  whether  TBRD  is 
calculated  above  or  below  the  horizon.   The  average  TBRD 
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Figure  21.  TBRD  Improvement- Starboard  Aspect 
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Figure  22.  TBRD  Improvement-Port  and  Stb  Aspect 
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Figure  23.  TBRD  Improvement-Bow  Aspect 
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improvemement  on  each  plot  is  calculated  only  for  the  below 
horizon  case  against  the  sea  water  background.   It  is 
expected  that  TBRD  against  the  sky  should  be  nearly  the  same 
in  both  polarizations.   This  data  does  not  contain  any 
meteorological,  target  temperature,  or  detailed  ship 
position  data  which  is  available  for  analysis  improvement. 

Very  general  trends  are  seen  for  TBRD  improvement,  but 
are  not  as  apparent  as  the  average  improvement  calculation 
numbers  in  the  plot.   Ship  size  varies  greatly  in  the  plots 
due  to  variation  in  the  numbers  of  pixels  recorded  in  the 
vertical  direction  and  is  most  pronounced  in  the  images 
taken  in  the  near  ranges.   It  is  seen  that  there  is  some 
polarization  effect  in  the  above  horizon  TBRD.   This  is 
probably  due  to  horizon  selecton  in  the  image  and  possible 
polarization  of  atmospheric  radiance  close  to  the  horizon. 
Analysis  was  not  performed  on  the  polarization  of  the  sky. 
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C.    FIELD  MEASUREMENTS  WITH  SPLIT-FIELD  POLARIMETER 

Preliminary  test  and  evaluation  imaging  was  done  on  3 
August,  1995  with  the  split-field  polarimeter  optical 
system.   The  target  is  the  USS  John  McCain,  an  Arleigh  Burke 
class  destroyer  present  in  Monterey  Bay  for  a  port  visit  at 
the  Naval  Postgraduate  School.   Thirteen  split-field  and  10 
test  images  were  taken  of  the  ship  and  black-body  calibra- 
tion sources.   The  target  was  at  a  stationary  range  of 
1.85km,  determined  by  target  geometrical  size.   All  images 
are  port  aspect.   All  of  the  images  taken  are  listed  by 
filename  in  Appendix  G.   Because  of  the  similarity  in  image 
composition,  only  one  image  will  be  analyzed  in  depth. 
Figure  24  shows  the  image  chosen  which  is  optimized  for 
focus,  level  gain  setting,  and  system  adjustment.    Images 
are  recorded  inverted  by  the  AGA  and  this  image  has  been 
artificially  flipped  upright  for  viewing.   The  top  image  is 
the  reflected  component  of  the  polarizing  splitter, 
horizontally  (s)  polarized,  and  the  bottom  is  the  trans- 
mitted vertically  (p)  polarized  image.   Consistent  with 
results  seen  for  the  internally  polarized  images  is  the 
increase  in  contrast  of  the  top  image.   Discussion  of 
specific  composition  of  each  channel  is  in  Chapter  VII. 
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Figure  25  shows  a  vertical  profile  taken  of  the  image 
through  the  hot  stack.   From  the  known  geometry  of  the  ship, 
no  significant  polarization  is  expected  for  the  stack  which 
is  just  aft  of  the  large  notch  in  the  top  of  the 
superstructure.   Polarization  of  emission  could  exist  in  the 
ship  panels  on  the  bow  as  they  slant  outward  toward  the 
horizontal  with  height.   The  bottom  of  the  profile 
corresponds  to  the  bottom  of  Fig.  24.   The  x-axis  denotes 


Figure  24.   NACIT  Split-field  polarized  image, 
USS  John  McCain,  3Aug95 
Top-Vertically  Polarized        TBRD=168 
Bottom-Horizontally  Polarized   TBRD=200 
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the  8  bit  digitized  value  of  radiant  intensity.   TBRD  for 
the  top  image  is  200  and  168  for  the  bottom.   This  gives  a 
TBRD  improvement  factor  =  1.19.   This  is  with  reference  to 
the  average  value  of  the  sea  water  background  taken  in  the 
foreground  along  the  azimuth  of  the  stack  in  the  middle. 
The  fog  above  the  image  is  assumed  unpolarized.   There  is  an 
overall  increase  in  image  intensity  in  the  vertically 
polarized  image  due  to  the  behavior  of  the  splitter  and  the 
background  plate  (Ch  7) . 
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Figure  25.   Vertical  Profile  of  Split  Image,  M0803011 
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The  added  reflectance  from  the  back  of  the  splitter  raises 
the  overall  level  of  the  vertically  polarized  image.   This 
is  another  justification  for  using  TBRD  to  compare  image 
contrast  differences. 

1.    Digital  Image  Subtraction 

The  goal  of  the  split  system  is  to  provide  simultaneous 
images  of  the  two  polarizations  which  allows  digital 
subtraction  of  the  images  to  visually  represent  TBRD 
improvement.   An  IDL  program  was  written  to  accomplish  a 
pixel  by  pixel  overlay  of  the  images  of  the  ships  (Appendix 
F) .    The  program  finds  the  'hottest'  pixel  in  the  image  and 
aligns  the  pixel  addresses  in  the  top  and  bottom  while 
subtracting  the  values.    The  resulting  image  displays  a 
value  of  0  (black)  where  the  images  had  the  same  intensity. 
The  areas  in  the  image  where  polarization  effects  occur  are 
the  brightest.   Figure  2  6  is  an  image  of  the  top  channel 
minus  the  bottom  channel  compensated  for  background 
contribution  to  the  bottom  image.   This  was  done  by 
adjusting  the  threshold  of  intensity  for  the  bottom  image. 
This  essentially  removes  any  polarization  effects  seen  on 
the  ship  because  of  their  small  effects. 

It  is  clearly  seen  that  the  water  shows  the  largest 
amount  of  difsference  in  the  image.   Irregularities  in 
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Figure  26.   Digitally  Subtracted  top  and  bottom 
split  image  channels.   M0803011,  3Aug95. 


in  the  image  plane  show  on  both  sides  of  the  subtracted 
image.   Technigues  are  under  development  to  balance  the 
intensity  of  the  two  images,  remove  background 
irregularites,  and  improve  alignment  of  the  two  images, 
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but  these  techniques  were  not  available  at  the  time 
these  images  were  processed.  The  images  are  shown  here  to 
illustrate  the  general  utility  and  feasibility. 
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IX.   DISCUSSION  AND  CONCLUSIONS 

The  split-field  polarizing  system  designed  and  con- 
structed in  this  work,  demonstrates  the  utility  of  simulta- 
neous  image  storing  for  processing,  while  preserving  image 
resolution. 

Subjective  and  numerical  analysis  of  the  data  from 
images  taken  with  the  split-field  system  demonstrated  excel- 
lent performance  in  image  quality,  contrast  improvement,  and 
utility.   Polarization  filtering  performance  for  a  sea  water 
background,  for  both  the  split-field  system  and  for  direct 
viewing  with  interchangeable  filters  in  the  AGEMA  imager, 
gave  results  comparable  to  those  obtained  in  the  MAPTIP 
experiments  where  an  external  BaFl  substrate  polarizer  was 
used.  This  was  determined  by  comparing  contrast  calculated 
in  terms  of  TBRD  from  previous  measurements  with  values  from 
the  split-field  and  internal  filtering  systems.   Both  the 
split-field  and  internally  filtered  direct  viewing  systems 
eliminated  the  narcissus  effects  that  degrade  the  use  of 
external  filters. 

Observations  of  the  target  scenes  show  quantitatively 
that: 

1.  Digital  subtraction  of  simultaneously   recorded  split- 
field  images  is  a  viable  technique  for  contrast  improvement. 
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2.  In  the  8-12  jam  region,  vertical  polarization  emission  of 
the  sea  surface  is  the  predominant  component  for  filtering. 

3.  Split-field  polarized  filtering  gives  preliminary  TBRD 
improvement  of  contrast  of  1.19.   This  is  driven  by  the 
splitter  performance  and  background  plate  effects. 

4.  Direct  viewing  with  internal  polarization  filtering  is 
capable  of  1.2-1.5  TBRD  improvement  of  contrast  for  sea- 
water  emission,  but  digital  subtraction  is  not  possible 
because  of  time  delay  between  acquisition  of  the 
orthogonally  polarized  images. 

5.  Both  the  split-field  and  direct  viewing  internal  filter- 
ing techniques  provide  contrast  improvement  on  the  order  of 
that  seen  in  the  MAPTIP  experiments. 

6.  The  combined  system  designed  and  constructed  for  this 
project  provides  dual  use  of  the  AGA-780  with  the  split- 
field  system  and  as  a  direct-view  internally-filtered  sys- 
tem. This  provides  two  highly  effective  polarizing  systems 
for  use  in  one  location  with  flexibility  in  system  selec- 
tion. 
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X.   RECOMMENDATIONS  FOR  FUTURE  WORK 

Further  work  needs  to  be  performed  in  all  aspects  of 
polarized  imaging  using  the  split  system  and  internal 
polarizing  systems.   In  March,  1996  an  Electro-Optical 
Propagation  Assessment  in  the  Coastal   Environment  (EOPACE) 
experiment  phase  will  take  place  near  San  Diego  CA.   In- 
depth  demonstration  and  field  implementation  of  both 
techniques  will  take  place  at  that  time.   Specific 
concentration  areas  for  both  systems  in  the  intermediate 
time  frame  include: 
Internal  Polarizers: 

-Quantify  the  spectral  transmissivity  experimentally 
for  each  polarizing  filter  to  verify  percent  polarization 
performance. 

-Re-calibrate  the  AGA  system  with  the  filters  installed 
prior  to  use  to  give  level  equalization  for  both  filters. 

Split-image  Polarimeter: 

-Conduct  laboratory  experiments  to  explore  and  remedy 
system  induced  gradient  across  the  f ield-of-view. 

-Realign  Split-field  optics  to  eliminate  any 
vignetting,  focussing  differences  between  the  two  channels, 
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and  edge  effects  due  to  optical  surface  misalignment. 

-Install  a  temperature  control  system  on  the  background 
plate  to  maintain  constant  radiance  contribution  to  the 
image  for  varying  environmental  conditions. 

-Insert  a  small  temperature  source  in  the  intermediate 
focal  plane  of  the  split-field  optics  to  give  a  constant 
calibration  source  for  real  time  split  channel  comparison. 

Overall  AGA-780  Thermal  Imaging  System: 

-Install  the  new  NEC  VERSA  P/75  75MHz  Pentium  Laptop 
computer  as  the  display  and  file  storage  server  to  replace 
the  current  CATS  38  6  computer. 

-Install  the  new  12  bit  A/D  converter  for  increasing 
system  dynamic  range  and  allowing  higher  data  resolution  for 
analysis. 
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APPENDIX  A. 
PERFORMANCE  CURVES  FOR  THE  POLARIZING  FILTERS 

This  appendix  contains  the  performance  curves  for  the 
polarizing  filters  used  in  the  internal  and  split  image 
configuration.   The  interpretation  of  the  data  is  included  in  the 
text  of  the  appropriate  chapters. 

The  data  were  provided  by  the  manufacturer  of  the 
polarizers,  Graseby-Specac.   The  measurements  were  taken  with  a 
spectrometer  using  the  actual  filter  shipped.   Where  two  filters 
were  required  for  the  crossed-grid  measurement  (trace  3) ,  another 
polarizer  nearly  identical  to  the  one  shipped  is  used. 

Trace  1-  transmittance  with  grid  perpendicular  to  the  E  vector  of 

incident  radiation. 

Trace  2-  transmittance  with  grid  parallel  to  the  E  vector  of 

incident  radiation. 

Trace  3-  transmittance  of  two  filters  with  grids  crossed. 

Vertical  Axis-  wavenumber/wavelength 
Horizontal  Axis-  transmittance 
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Figure  Al .   Performance  curves  for  50mm  Ge  Splitter 
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Figure  A2 .   Performance  curves  for  KRS-5  9.5mm  polarizer  (#75) 
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Figure  A3.   Performance  curves  for  KRS-5  9.5mm  polarizer  (#74) 


90 


APPENDIX  B. 
SPECTRAL  TRANSMITTANCE /REFLECTANCE  OF  THE  50  MM  POLARIZER 


This  appendix  contains  the  Perkin-Elmer 
spectophotometer  traces  for  the  50mm  polarizing  splitter. 
In  Fig.  Bl,  the  traces  with  the  x  designation  are  for  the 
polarizer  oriented  45  degrees  to  the  incoming  beam.   The  K1N 
trace  is  the  combination  of  the  traces  and  gives  the  overall 
transmittance  of  the  splitter  normal  to  the  incoming  beam. 
The  x  symbols  along  the  K1N  trace  are  the  overall 
transmittance  of  the  splitter  at  45  degrees  to  the  incoming 
beam.   It  is  seen  that  the  peak  occurs  at  a  slightly  shorter 
wavelength  than  at  normal  incidence. 

In  Fig.B2,  K1R  is  the  overall  reflection  efficiency  in 
the  configuration  described  in  Chapt .  VI. 
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Figure  Bl .   Polarization  Transmittance  of  50mm  Polarizer 
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Figure  B2 .   Reflection  of  Polarized  Radiation  by  the  Grid. 
Transmittance  axis  denotes  Reflectance  (%)  . 
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APPENDIX  C. 

MEASURED  RADIANCE  VALUES  FOR  THE  NPS-NACIT  SPLIT-FIELD 

POLARIMETER 


The  following  data  were  taken  as  input  to  the  radiance 
analysis  of  polarizing  splitter  performance. 

Lab  Conditions: 

Room  Temp  17C, 

Source-Agema  Black-Body  @60C 

Polarized  Source-BaFl  IR  Polarizer  Pol%~99% 

Sensor-AGA  780  LW  Band,  f/1.87,  Range  2 

Background  Plate  Temp-  ~0C  (heavy  frost  surface) 

£  =21.7   g  =23.2 


=4.2   g  =11.2 
=8.6   £_  =5.4 


I  =34.7  I   =51.0 
1^=50.3   1^=36.9 
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APPENDIX  D. 
CALIBRATION  CURVES  FOR  FILTERED  AND  UNFILTERED  AGA-780 
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Figure  D-l. 
present . 


Calibration  for  the  AGA  with  no  filters 
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Figure  D2 .   Calibration  Data  for  AGA  with  internal  filter. 
Pos  5  denotes  one  of  the  small  internal  filters  located  in 
position  5  of  the  AGA  filter  wheel.   Position  4  is  the  other 
with  orthogonal  filter  orientation.   A  small  difference  is 
seen  in  thermal  level  with  temperature  between  the  two 
filters . 
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Figure  D3 .   Calibration  curve  for  the  Split-field 
Polarimeter .  This  shows  level  reading  with  temperature  for 
the  two  channels  of  the  split-field  polarimeter.   LWPOL  is 
for  the  large  BaFl  external  polarizer  for  comparison.   Upper 
spot  denotes  the  splitter  transmitted  channel  imaged  by  the 
AGA. 
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APPENDIX  E. 

COMPILATION  OF  MOSS  LANDING  IMAGES 

MAY,  1995 


The  following  is  a  complete  image  compilation  of  the  entire 
Moss  Landing  experiment.   The  data  tabulated  is  incomplete  and 
only  shows  what  is  necessary  for  preliminary  TBRD  comparison  of 
images.   A  sheet  depicting  a  number  of  representitive  images  is 
appended  to  this  list.   Images  corresponding  to  each  of  the  table 
entries  are  available  in  digital  form  on  hard  disk  at  the  NACIT 
laboratory. 
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The  title  header  gives  data  about  each  of  the  images  taken. 
Pol-polarization  orientation  of  the  image,  h, v-horiz, vert .  x-image  not 

analyzable  for  one  reason  (see  image).   u-unpolarized  image. 
Level- thermal  level  setting  of  the  AGA  control  knob. 
Range-thermal  level  range  setting  of  the  AGA  control. 

Offset-Offset  of  the  calibration  constant  for  image  temperature  compensation. 
Ovrflw,  Undrflw-number  of  pixels  in  the  image  with  radiance  values  above  or 
below  the  0-255  range. 
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level 
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M5L19092.IMG  5/19/95 
M5L19093.IMG  5/19/95 
M5L19094.IMG  5/19/95 
M5L19095.IMG  5/19/95 
M5L19096.IMG  5/19/95 
M5L19097.IMG  5/19/95 
M5L19098.IMG  5/19/95 
M5L19099.IMG  5/19/95 
M5L19100.IMG  5/19/95 
M5L19101.IMG  5/19/95 
M5L19102.IMG  5/19/95 
M5L19103.IMG  5/19/95 
M5L19104.IMG  5/19/95 
M5L19105.IMG  5/19/95 
M5L19106.IMG  5/19/95 
M5L19107.IMG  5/19/95 
M5L19108.IMG  5/19/95 
M5L19109.IMG  5/19/95 
M5L19110.IMG  5/19/95 
M5L19111.IMG  5/19/95 
M5L19112.IMG  5/19/95 
M5L19113.IMG  5/19/95 
M5L19114.IMG  5/19/95 
M5L19115.IMG  5/19/95 
M5L19116.IMG  5/19/95 
M5L19117.IMG  5/19/95 
M5L19118.IMG  5/19/95 
M5L19119.IMG  5/19/95 
M5L19120.IMG  5/19/95 
M5L19121.IMG  5/19/95 
M5L19122.IMG  5/19/95 
M5L19123.IMG  5/19/95 
M5L19124.IMG  5/19/95 
M5L19125.IMG  5/19/95 
M5L19126.IMG  5/19/95 
M5L19127.IMG  5/19/95 
M5L19128.IMG  5/19/95 
M5L19129.IMG  5/19/95 
M5L19130.IMG  5/19/95 
M5L19131.IMG  5/19/95 
M5L19132.IMG  5/19/95 
M5L19133.IMG  5/19/95 
M5L19134.IMG  5/19/95 
M5L19135.IMG  5/19/95 
M5L19136.IMG  5/19/95 
M5L19137.IMG  5/19/95 
M5L19138.IMG  5/19/95  10:04:15 
M5L19139.IMG  5/19/95  10:17:50 
M5L19140.IMG  5/19/95  10:17:54 
M5L19141.IMG  5/19/95  10:18:03 
M5L19142.IMG  5/19/95  10:18:30 
M5L19143.IMG  5/19/95  10:18:59 
M5L19144.IMG  5/19/95  10:19:04 
M5L19145.IMG  5/19/95  10:19:14 
M5L19146.IMG  5/19/95  10:19:21 
M5L19147.IMG  5/19/95  10:35:59 
M5L19148.IMG  5/19/95  10:37:23 
M5L19149.IMG  5/19/95  10:37:50 
M5L19150.IMG  5/19/95  10:39:50 
M5L19151.IMG  5/19/95  10:39:55 
M5L19152.IMG  5/19/95  10:40:04 
M5L19153.IMG  5/19/95  10:41:01 
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55:04 
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44.00 
47.10 
47.10 
44.20 
48.60 
48.10 
43.70 
46.90 
46.90 
43.50 
47.10 
47.60 
43.50 
47.60 
47.60 
44.00 
47.90 
47.90 
47.90 
45.20 
47.40 
47.40 
43.70 
48.10 
48.10 
44.40 
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44.00 
47.90 
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47.10 
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M5L19154.IMG  5/19/95  10:41:13  v  48.10  5 

M5L19155.IMG  5/19/95  10:41:24  u  44.20  5, 

M5L19156.IMG  5/19/95  10:44:33  h  46.60  2, 

M5L19157.IMG  5/19/95  10:44:36  v  46.60  2, 

M5L19158.IMG  5/19/95  10:44:46  u  44.00  2, 

M5L19159.IMG  5/19/95  10:48:01  h  47.60  2. 

M5L19160.IMG  5/19/95  10:48:06  v  47.60  2. 

M5L19161.IMG  5/19/95  10:48:20  u  46.20  2. 

M5L19162.IMG  5/19/95  10:49:47  h  46.60  2, 

M5L19163.IMG  5/19/95  10:49:53  v  46.60  2. 

M5L19164.IMG  5/19/95  10:50:01  u  43.20  2. 

M5L19165.IMG  5/19/95  10:51:40  h  46.90  2. 

M5L19166.IMG  5/19/95  10:51:43  v  46.90  2. 

M5L19167.IMG  5/19/95  10:51:49  u  43.50  2. 

M5L19168.IMG  5/19/95  10:52:09  h  47.60  5. 

M5L19169.IMG  5/19/95  10:52:14  v  47.60  5. 

M5L19170.IMG  5/19/95  10:52:22  u  44.00  5. 

M5L19171.IMG  5/19/95  10:53:53  h  46.60  2. 

M5L19172.IMG  5/19/95  10:53:56  v  46.60  2. 

M5L19173.IMG  5/19/95  10:54:03  u  43.50  2. 

M5L19174.IMG  5/19/95  10:57:32  h  46.20  2. 

M5L19175.IMG  5/19/95  10:57:36  v  46.20  2. 

M5L19176.IMG  5/19/95  10:57:46  u  43.70  2. 

M5L19177.IMG  5/19/95  10:58:02  h  47.40  5. 

M5L19178.IMG  5/19/95  10:58:06  v  47.40  5. 

M5L19179.IMG  5/19/95  10:58:13  u  44.00  5. 

M5L19180.IMG  5/19/95  10:59:45  h  46.60  2. 

M5 LI 9181. IMG  5/19/95  10:59:49  v  46.60  2. 

M5L19182.IMG  5/19/95  10:59:57  u  43.20  2. 

M5L19183.IMG  5/19/95  11:00:20  h  46.20  5. 

M5L19184.IMG  5/19/95  11:00:23  v  46.20  5. 

M5L19185.IMG  5/19/95  11:00:30  u  44.20  5. 

M5L19186.IMG  5/19/95  11:01:46  h  46.90  2. 

M5L19187.IMG  5/19/95  11:01:49  v  46.90  2. 

M5L19188.IMG  5/19/95  11:01:56  u  43.70  2. 

M5L19189.IMG  5/19/95  11:02:09  h  45.90  5. 

M5L19190.IMG  5/19/95  11:02:14  v  45.90  5. 

M5L19191.IMG  5/19/95  11:02:20  u  44.70  5. 

M5L19192.IMG  5/19/95  11:04:24  h  46.60  2. 

M5L19193.IMG  5/19/95  11:04:28  v  46.60  2. 

M5L19194.IMG  5/19/95  11:04:39  u  42.70  2. 

M5L19195.IMG  5/19/95  11:05:06  h  47.10  5. 

M5L19196.IMG  5/19/95  11:05:09  v  47.10  5. 

M5L19197.IMG  5/19/95  11:05:16  u  44.70  5. 

M5L19198.IMG  5/19/95  11:06:01  h  48.40  10. 

M5L19199.IMG  5/19/95  11:06:04  v  48.40  10. 

M5L19200.IMG  5/19/95  11:06:10  u  45.90  10. 

M5L19201.IMG  5/19/95  11:08:23  h  46.90  2. 

M5L19202.IMG  5/19/95  11:08:26  v  46.90  2. 

M5L19203.IMG  5/19/95  11:08:35  u  43.20  2. 

M5L19204.IMG  5/19/95  11:08:49  h  45.70  5. 

M5L19205.IMG  5/19/95  11:08:57  v  46.40  5. 

M5L19206.IMG  5/19/95  11:09:04  u  44.70  5. 

M5L19207.IMG  5/19/95  11:09:15  h  46.20  10. 

M5L19208.IMG  5/19/95  11:09:19  v  46.20  10. 

M5L19209.IMG  5/19/95  11:09:25  u  46.60  10. 

M5L19210.IMG  5/19/95  11:11:11  h  48.80  5. 

M5L19211.IMG  5/19/95  11:11:17  v  48.80  5. 

M5L19212.IMG  5/19/95  11:11:28  u  46.60  5. 

M5L19213.IMG  5/19/95  11:11:43  u  46.60  5. 

M5L19214.IMG  5/19/95  11:14:36  h  46.90  2. 

M5L19215.IMG  5/19/95  11:14:39  v  46.90  2. 


-5. 

0. 

0. 

-5. 

0. 

0. 

-5. 

0. 

0. 

-5. 

0. 

0. 

-5. 

0. 

0. 

-5. 

19. 

0. 

-5. 

142. 

0. 

■5. 

334. 

1 

■5. 

0. 

0. 

-5. 

0. 

0. 

■5. 

1. 

0. 

■5. 

0. 

0. 

-5. 

0. 

0. 

-5. 

0. 

0. 

-5. 

0. 

0. 

•5. 

0. 

0. 

•5. 

0. 

0. 

-5. 

0. 

0. 

■5. 

0. 

0. 

■5. 

0. 

0. 

■5. 

0. 

0. 

■5. 

0. 

0. 

■5. 

0. 

0. 

■5. 

0. 

0. 

■5. 

0. 

0. 

■5. 

0. 

0. 

■5. 

0. 

0. 

■5. 

0. 

0. 

■5. 

71. 

0. 

•5. 

0. 

0. 

■5. 

0. 

0. 

■5. 

0. 

0. 

■5. 

2. 

0. 

■5. 

21. 

0. 

■5. 

111. 

0. 

5. 

0. 

0. 

5. 

0. 

0. 

5. 

11. 

0. 

5. 

28. 

0. 

5. 

41. 

0. 

5. 

768. 

0. 

5. 

0. 

0. 

5. 

0. 

0. 

5. 

24. 

0. 

5. 

0. 

0. 

5. 

0. 

0. 

5. 

0. 

0. 

5. 

2. 

0. 

5. 

3. 

0. 

5. 

186. 

0. 

5. 

0. 

0. 

5. 

0. 

0. 

5. 

0. 

0. 

5. 

0. 

0. 

5. 

0. 

0. 

5. 

0. 

0. 

5. 

12. 

3681 

5. 

15. 

802. 

5. 

338. 

5606 

5. 

1. 

5640 

5. 

0. 

0. 

5. 

0. 

0. 

13059, 


108 


M5L19216.IMG  5/19/95  11:14:46  u  43.50  2 

M5L19217.IMG  5/19/95  11:15:05  h  47.10  5 

M5L19218.IMG  5/19/95  11:15:08  v  47.10  5 

M5L19219.IMG  5/19/95  11:15:14  u  44.90  5, 

M5L19220.IMG  5/19/95  11:15:39  h  48.10  10, 

M5L19221.IMG  5/19/95  11:15:42  v  48.10  10. 

M5L19222.IMG  5/19/95  11:15:49  u  46.40  10. 

M5L19223.IMG  5/19/95  11:16:53  h  46.90  2. 

M5L19224.IMG  5/19/95  11:16:56  v  46.90  2. 

M5L19225.IMG  5/19/95  11:17:01  u  43.50  2. 

M5L19226.IMG  5/19/95  11:17:27  h  46.90  5. 

M5L19227.IMG  5/19/95  11:17:30  v  46.90  5. 

M5L19228.IMG  5/19/95  11:17:35  u  44.70  5. 

M5L19229.IMG  5/19/95  11:17:53  h  48.40  10. 

M5L19230.IMG  5/19/95  11:17:56  v  48.40  10. 

M5L19231.IMG  5/19/95  11:18:01  u  46.20  10. 

M5L19232.IMG  5/19/95  11:19:11  h  47.40  2. 

M5L19233.IMG  5/19/95  11:19:14  v  46.90  2. 

M5L19234.IMG  5/19/95  11:19:18  u  43.70  2. 

M5L19235.IMG  5/19/95  11:19:38  h  48.80  5. 

M5L19236.IMG  5/19/95  11:19:41  v  48.80  5. 

M5L19237.IMG  5/19/95  11:19:48  u  44.40  5. 

M5L19238.IMG  5/19/95  11:20:07  h  48.80  10. 

M5L19239.IMG  5/19/95  11:20:10  v  48.80  10. 

M5L19240.IMG  5/19/95  11:20:19  u  46.20  10. 

M5L19241.IMG  5/19/95  11:21:23  h  46.90  2. 

M5L19242.IMG  5/19/95  11:21:26  v  46.90  2. 

M5L19243.IMG  5/19/95  11:21:29  u  43.20  2. 

M5L19244.IMG  5/19/95  11:21:48  h  48.40  5. 

M5L19245.IMG  5/19/95  11:21:51  v  48.40  5. 

M5L19246.IMG  5/19/95  11:21:56  u  44.90  5. 

M5L19247.IMG  5/19/95  11:22:16  h  48.60  10. 

M5L19248.IMG  5/19/95  11:22:19  v  48.60  10. 

M5L19249.IMG  5/19/95  11:22:26  u  46.20  10. 

M5L19250.IMG  5/19/95  11:22:58  h  50.80  20. 

M5L19251.IMG  5/19/95  11:23:01  v  50.80  20. 

M5L19252.IMG  5/19/95  11:23:08  u  50.30  20. 

M5L19253.IMG  5/19/95  11:24:50  u  44.90  5. 

M5L19254.IMG  5/19/95  11:25:29  h  48.60  5. 

M5L19255.IMG  5/19/95  11:25:32  v  48.60  5. 

M5L19256.IMG  5/19/95  11:25:36  u  44.70  5. 

M5L19257.IMG  5/19/95  11:25:50  h  49.60  10. 

M5L19258.IMG  5/19/95  11:25:53  v  49.60  10. 

M5L19259.IMG  5/19/95  11:25:58  u  46.60  10. 

M5L19260.IMG  5/19/95  11:26:14  h  47.60  20. 

M5L19261.IMG  5/19/95  11:26:17  v  47.60  20. 

M5L19262.IMG  5/19/95  11:26:20  u  47.60  20. 

M5L19263.IMG  5/19/95  11:28:22  h  48.80  5. 

M5L19264.IMG  5/19/95  11:28:25  v  48.80  5. 

M5L19265.IMG  5/19/95  11:28:33  u  44.90  5. 

M5L19266.IMG  5/19/95  11:28:56  h  46.90  10. 

M5L19267.IMG  5/19/95  11:28:59  v  46.90  10. 

M5L19268.IMG  5/19/95  11:29:02  u  46.90  10. 

M5L19269.IMG  5/19/95  11:30:11  h  47.90  5. 

M5L19270.IMG  5/19/95  11:30:14  v  47.90  5. 

M5L19271.IMG  5/19/95  11:30:24  u  44.20  5. 

M5L19272.IMG  5/19/95  11:30:40  h  48.10  10. 

M5L19273.IMG  5/19/95  11:30:44  v  48.10  10. 

M5L19274.IMG  5/19/95  11:30:51  u  47.40  10. 

M5L19275.IMG  5/19/95  11:32:15  h  47.90  5. 

M5L19276.IMG  5/19/95  11:32:18  v  47.90  5. 

M5L19277.IMG  5/19/95  11:32:22  u  44.90  5. 
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M5L19278.IMG  5/19/95  11:32:38 
M5L19279.IMG  5/19/95  11:32:41 
M5L19280.IMG  5/19/95  11:32:44 
M5L19281.IMG  5/19/95  11:37:53 
M5L19282.IMG  5/19/95  11:37:58 
M5L19283.IMG  5/19/95  11:38:05 
M5L19284.IMG  5/19/95  11:38:32 
M5L19285.IMG  5/19/95  11:38:35 
M5L19286.IMG  5/19/95  11:38:40 
M5L19287.IMG  5/19/95  11:40:25 
M5L19288.IMG  5/19/95  11:40:50 
M5L19289.IMG  5/19/95  11:40:59 
ML516001.IMG  5/15/95  11:48:25 
ML516002.IMG  5/15/95  11:48:48 
ML516010.IMG  5/15/95  11:49:49 
ML516011.IMG  5/15/95  11:50:29 
ML516012.IMG  5/15/95  11:50:56 
ML516013.IMG  5/15/95  11:50:56 
ML516014.IMG  5/15/95  11:50:56 
ML516015.IMG  5/15/95  12:02:30 
ML516016.IMG  5/16/95  9:29:08 
ML516017.IMG  5/16/95 
ML516018.IMG  5/16/95 
ML516019.IMG  5/16/95 
ML516020.IMG  5/16/95 
ML516021.IMG  5/16/95 
ML516022.IMG  5/16/95 
ML516023.IMG  5/16/95 
ML516024.IMG  5/16/95 
ML516025.IMG  5/16/95 
ML516026.IMG  5/16/95 
ML516027.IMG  5/16/95 
ML516028.IMG  5/16/95 
ML516029.IMG  5/16/95 
ML516030.IMG  5/16/95 
ML516031.IMG  5/16/95 
ML516032.IMG  5/16/95 
ML516033.IMG  5/16/95  10:00:57 
ML516034.IMG  5/16/95  10:01:10 
ML516035.IMG  5/16/95  10:01:30 
ML516036.IMG  5/16/95  10:01:56 
ML516037.IMG  5/16/95  10:03:05 
ML516038.IMG  5/16/95  10:03:16 
ML516039.IMG  5/16/95  10:03:42 
ML516040.IMG  5/16/95  10:04:23 
ML516041.IMG  5/16/95  10:04:35 
ML516042.IMG  5/16/95  10:04:56 
ML516043.IMG  5/16/95  10:05:09 
ML516044.IMG  5/16/95  10:06:26 
ML516045.IMG  5/16/95  10:06:35 
ML516046.IMG  5/16/95  10:07:09 
ML516047.IMG  5/16/95  10:08:09 
ML516048.IMG  5/16/95  10:08:15 
ML516049.IMG  5/16/95  10:08:26 
ML516050.IMG  5/16/95  10:09:01 
ML516051.IMG  5/16/95  10:09:52 
ML516052.IMG  5/16/95  10:10:02 
ML516053.IMG  5/16/95  10:10:26 
ML516054.IMG  5/16/95  10:10:33 
ML516055.IMG  5/16/95  10:11:17 
ML516056.IMG  5/16/95  10:13:02 
ML516057.IMG  5/16/95  10:13:12 
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ML516058.IMG  5/16/95  10:13:28 
ML516059.IMG  5/16/95  10:15:43 
ML516060.IMG  5/16/95  10:15:56 
ML516061.IMG  5/16/95  10:16:26 
ML516062.IMG  5/16/95  10:16:41 
ML516063.IMG  5/16/95  10:18:43 
ML516064.IMG  5/16/95  10:18:52 
ML516065.IMG  5/16/95  10:19:15 
ML516066.IMG  5/16/95  10:19:58 
ML516067.IMG  5/16/95  10:20:06 
ML516068.IMG  5/16/95  10:20:16 
ML516069.IMG  5/16/95  10:20:54 
ML516070.IMG  5/16/95  10:21:00 
ML516071.IMG  5/16/95  10:21:11 
ML516072.IMG  5/16/95  10:23:14 
ML516073.IMG  5/16/95  10:23:23 
ML516074.IMG  5/16/95  10:23:33 
ML516075.IMG  5/16/95  10:24:16 
ML516076.IMG  5/16/95  10:24:26 
ML516077.IMG  5/16/95  10:24:36 
ML516078.IMG  5/16/95  10:25:15 
ML516079.IMG  5/16/95  10:25:22 
ML516080.IMG  5/16/95  10:25:30 
ML516081.IMG  5/16/95  10:27:54 
ML516082.IMG  5/16/95  10:28:32 
ML516083.IMG  5/16/95  10:28:56 
ML516084.IMG  5/16/95  10:29:21 
ML516085.IMG  5/16/95  10:30:11 
ML516086.IMG  5/16/95  10:30:35 
ML516087.IMG  5/16/95  10:32:59 
ML516088.IMG  5/16/95  10:33:13 
ML516089.IMG  5/16/95  10:33:36 
ML516090.IMG  5/16/95  10:33:52 
ML516091.IMG  5/16/95  10:35:03 
ML516092.IMG  5/16/95  10:35:16 
ML516093.IMG  5/16/95  10:35:40 
ML516094.IMG  5/16/95  10:35:53 
ML516095.IMG  5/16/95  10:37:31 
ML516096.IMG  5/16/95  10:37:46 
ML516097.IMG  5/16/95  10:37:59 
ML516098.IMG  5/16/95  10:38:29 
ML516099.IMG  5/16/95  10:38:38 
ML516100.IMG  5/16/95  10:38:55 
ML516101.IMG  5/16/95  10:39:13 
ML516102.IMG  5/16/95  10:41:36 
ML516103.IMG  5/16/95  10:41:49 
ML516104.IMG  5/16/95  10:42:19 
ML516105.IMG  5/16/95  10:42:34 
ML516106.IMG  5/16/95  10:43:12 
ML516107.IMG  5/16/95  10:43:22 
ML516108.IMG  5/16/95  10:43:42 
ML516109.IMG  5/16/95  10:45:07 
ML516110.IMG  5/16/95  10:45:16 
ML516111.IMG  5/16/95  10:45:37 
ML516112.IMG  5/16/95  10:47:01 
ML516113.IMG  5/16/95  10:47:32 
ML516114.IMG  5/16/95  10:48:32 
ML516115.IMG  5/16/95  10:48:41 
ML516116.IMG  5/16/95  10:48:52 
ML516117.IMG  5/16/95  10:50:30 
ML516118.IMG  5/16/95  10:50:44 
ML516119.IMG  5/16/95  10:51:02 
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X 

0 

00 

2 

ML522003 

IMG 

5/22/95 

8 

12 

45 

X 

0 

00 

2 

ML522004 

IMG 

5/22/95 

8 

12 

53 

X 

0 

00 

2 

ML522005 

IMG 

5/22/95 

8 

46 

39 

X 

44 

00 

5 

ML522006 

IMG 

5/22/95 

8 

47 

30 

u 

43 

70 

2 

ML522007 

IMG 

5/22/95 

8 

48' 

05 

h 

45 

90 

2 

ML522008 

IMG 

5/22/95 

8 

48 

20 

V 

46 

20 

2 

ML522009 

IMG 

5/22/95 

8 

49 

06 

u 

44 

90 

5 

ML522010 

IMG 

5/22/95 

8 

49: 

51 

h 

46 

20 

2 

ML522011 

IMG 

5/22/95 

8 

50: 

07 

V 

46 

20 

2 

ML522012 

IMG 

5/22/95 

8 

50. 

46 

u 

44 

40 

5 

ML522013 

IMG 

5/22/95 

8 

52: 

25 

h 

45 

70 

2 

ML522014 

IMG 

5/22/95 

8 

52: 

47 

V 

46 

20 

2 

ML522015 

IMG 

5/22/95 

8 

53: 

00 

u 

44 

90 

5 

ML522016 

IMG 

5/22/95 

8 

59- 

23 

h 

45 

70 

2 

ML522017. 

IMG 

5/22/95 

8 

59: 

30 

V 

45 

70 

2 

ML522018. 

IMG 

5/22/95 

8 

59: 

42 

u 

44 

90 

5 

ML522019 

IMG 

5/22/95 

9 

00: 

06 

h 

45 

90 

2 

ML522020. 

IMG 

5/22/95 

9 

00: 

11 

V 

45 

90 

2 

ML522021. 

IMG 

5/22/95 

9 

00: 

21 

u 

44 

90 

5 

ML522022 

IMG 

5/22/95 

9 

01: 

38 

u 

44 

20 

5 

ML522023 

IMG 

5/22/95 

9 

01: 

58 

h 

46 

40 

2 

ML522024 

IMG 

5/22/95 

9 

02: 

03 

V 

46 

40 

2 

ML522025 

IMG 

5/22/95 

9 

02: 

11 

u 

45 

40 

5 

ML522026 

IMG 

5/22/95 

9 

03: 

44 

h 

45 

70 

2 

ML522027 

IMG 

5/22/95 

9 

03: 

50 

V 

45 

70 

2 

ML522028 

IMG 

5/22/95 

9 

03: 

59 

u 

44 

00 

5 

ML522029 

IMG 

5/22/95 

9 

04: 

21 

h 

45 

70 

2 

ML522030 

IMG 

5/22/95 

9 

04. 

25 

V 

45. 

70 

2 

ML522031 

IMG 

5/22/95 

9 

04: 

39 

u 

44 

00 

5 

ML522032 

IMG 

5/22/95 

9 

05. 

22 

h 

45. 

20 

5 

ML522033 

IMG 

5/22/95 

9 

05: 

26 

V 

45 

20 

5 

ML522034 

IMG 

5/22/95 

9 

05: 

32 

u 

45. 

20 

5 

0. 

0. 

19081. 

0. 

195. 

0. 

0. 

327. 

0. 

0. 

253. 

0. 

0. 

0. 

19040. 

0. 

163. 

0. 

0. 

285. 

0. 

0. 

160. 

0. 

0. 

0. 

19082. 

0. 

213. 

0. 

0. 

350. 

0. 

0. 

221. 

0. 

0. 

0. 

19097. 

0. 

155. 

0. 

0. 

282. 

0. 

0. 

308. 

0. 

0. 

276. 

0. 

0. 

0. 

19162. 

0. 

996. 

0. 

0. 

545. 

0. 

0. 

323. 

0. 

0. 

0. 

19088. 

0. 

213. 

0. 

0. 

0. 

18994. 

0. 

316. 

0. 

0. 

0. 

19077. 

0. 

590. 

1481 

0. 

19600.      0 

0. 

19600.      0 

0. 

19600.      0 

0. 

19600.      0 

3. 

57. 

0. 

5. 
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5. 
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0. 

5. 
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0. 

5. 
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5. 

40. 

2. 

5. 
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5. 
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5. 
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5. 
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5. 
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0. 

5. 

26. 
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5. 
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5. 

5. 
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0. 

5. 
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5. 
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5. 
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5. 
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5. 
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9117. 

5. 

0. 
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0. 
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5. 

0. 
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5. 

0. 

0. 

5. 

0. 

5. 

5. 

0. 

0. 

5. 

0. 

0. 

5. 

0. 

3312. 
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ML522035.IMG  5/22/95  9:07:15  h  45.40  2 

ML522036.IMG  5/22/95  9:07:22  v  45.40  2, 

ML522037.IMG  5/22/95  9:07:31  u  44.40  5. 

ML522038.IMG  5/22/95  9:07:52  h  45.70  2, 

ML522039.IMG  5/22/95  9:07:56  v  45.70  2. 

ML522040.IMG  5/22/95  9:08:06  u  44.00  5. 

ML522041.IMG  5/22/95  9:08:32  h  45.90  5. 

ML522042.IMG  5/22/95  9:08:35  v  45.90  5. 

ML522043.IMG  5/22/95  9:08:47  u  44.20  5. 

ML522044.IMG  5/22/95  9:16:02  h  45.90  2. 

ML522045.IMG  5/22/95  9:16:06  v  45.90  2. 

ML522046.IMG  5/22/95  9:16:18  u  44.40  5. 

ML522047.IMG  5/22/95  9:16:40  h  45.70  2. 

ML522048.IMG  5/22/95  9:16:44  v  45.70  2. 

ML522049.IMG  5/22/95  9:16:53  u  44.90  5. 

ML522050.IMG  5/22/95  9:19:01  h  45.70  2. 

ML522051.IMG  5/22/95  9:19:06  v  45.70  2. 

ML522052.IMG  5/22/95  9:19:15  u  44.00  5. 

ML522053.IMG  5/22/95  9:19:31  h  45.90  2. 

ML522054.IMG  5/22/95  9:19:34  v  45.90  2. 

ML522055.IMG  5/22/95  9:19:41  u  44.20  5. 

ML522056.IMG  5/22/95  9:20:49  h  45.40  5. 

ML522057.IMG  5/22/95  9:21:01  v  45.40  5. 

ML522058.IMG  5/22/95  9:21:15  u  44.20  5. 

ML522059.IMG  5/22/95  9:22:23  h  45.40  2. 

ML522060.IMG  5/22/95  9:22:29  v  45.40  2. 

ML522061.IMG  5/22/95  9:22:39  u  44.40  5. 

ML522062.IMG  5/22/95  9:22:51  h  45.90  2. 

ML522063.IMG  5/22/95  9:22:55  v  45.90  2. 

ML522064.IMG  5/22/95  9:23:02  u  44.20  5. 

ML522065.IMG  5/22/95  9:26:48  h  45.90  2. 

ML522066.IMG  5/22/95  9:26:52  v  45.90  2. 

ML522067.IMG  5/22/95  9:29:34  h  46.40  2. 

ML522068.IMG  5/22/95  9:29:40  v  46.40  2. 

ML522069.IMG  5/22/95  9:29:51  u  44.90  5. 

ML522070.IMG  5/22/95  9:30:35  h  46.20  2. 

ML522071.IMG  5/22/95  9:30:39  v  46.20  2. 

ML522072.IMG  5/22/95  9:30:49  u  44.70  5. 

ML522073.IMG  5/22/95  9:32:37  h  46.20  2. 

ML522074.IMG  5/22/95  9:32:41  v  46.20  2. 

ML522075.IMG  5/22/95  9:32:49  u  44.40  5. 

ML522076.IMG  5/22/95  9:33:05  h  46.40  2. 

ML522077.IMG  5/22/95  9:33:11  v  46.40  2. 

ML522078.IMG  5/22/95  9:33:23  u  44.20  5. 

ML522079.IMG  5/22/95  9:34:45  h  45.70  2. 

ML522080.IMG  5/22/95  9:34:49  v  45.70  2. 

ML522081.IMG  5/22/95  9:34:56  u  44.20  5. 

ML522082.IMG  5/22/95  9:35:11  h  46.40  2. 

ML522083.IMG  5/22/95  9:35:14  v  46.40  2. 

ML522084.IMG  5/22/95  9:35:21  u  44.20  5. 

ML522085.IMG  5/22/95  9:45:03  h  46.40  2. 

ML522086.IMG  5/22/95  9:45:07  v  46.40  2. 

ML522087.IMG  5/22/95  9:45:15  u  44.90  5. 

ML522088.IMG  5/22/95  9:45:39  h  46.60  2. 

ML522089.IMG  5/22/95  9:45:43  v  46.60  2. 

ML522090.IMG  5/22/95  9:45:53  u  43.70  5. 

ML522091.IMG  5/22/95  9:47:46  h  46.20  2. 

ML522092.IMG  5/22/95  9:48:04  h  46.20  2. 

ML522093.IMG  5/22/95  9:48:12  v  46.20  2. 

ML522094.IMG  5/22/95  9:48:22  u  44.00  5. 

ML522095.IMG  5/22/95  9:48:46  x  46.40  2. 

ML522096.IMG  5/22/95  9:48:51  h  46.40  2. 
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ML522097.IMG  5/22/95  9:49:27  v  46.20  2 

ML522098.IMG  5/22/95  9:49:31  u  46.20  2 

ML522099.IMG  5/22/95  9:49:45  u  44.00  2 

ML522100.IMG  5/22/95  9:51:20  h  45.90  2 

ML522101.IMG  5/22/95  9:51:25  v  45.90  2 

ML522102.IMG  5/22/95  9:51:37  u  43.70  2 

ML522103.IMG  5/22/95  9:51:53  h  46.20  2 

ML522104.IMG  5/22/95  9:51:57  v  46.20  2 

ML522105.IMG  5/22/95  9:52:07  u  44.20  2 

ML522106.IMG  5/22/95  10:36:07  h  47.10  2, 

ML522107.IMG  5/22/95  10:36:15  v  47.10  2 

ML522108.IMG  5/22/95  10:36:32  v  46.90  2. 

ML522109.IMG  5/22/95  10:36:51  u  44.40  2 

ML522110.IMG  5/22/95  10:37:20  h  46.90  2. 

ML522111.IMG  5/22/95  10:37:25  v  46.40  2. 

ML522112.IMG  5/22/95  10:37:35  u  44.70  2. 

ML522113.IMG  5/22/95  10:38:29  h  46.90  2. 

ML522114.IMG  5/22/95  10:38:33  v  46.90  2. 

ML522115.IMG  5/22/95  10:38:52  u  44.20  2. 

ML522116.IMG  5/22/95  10:39:25  h  46.90  2. 

ML522117.IMG  5/22/95  10:39:29  v  46.90  2. 

ML522118.IMG  5/22/95  10:39:41  u  44.20  2. 

ML522119.IMG  5/22/95  10:41:09  h  46.90  2, 

ML522120.IMG  5/22/95  10:41:14  v  46.90  2. 

ML522121.IMG  5/22/95  10:41:25  u  44.20  2, 

ML522122.IMG  5/22/95  10:41:40  h  47.10  2. 

ML522123.IMG  5/22/95  10:41:45  v  47.10  2. 

ML522124.IMG  5/22/95  10:41:56  u  44.40  2. 

ML522125.IMG  5/22/95  10:50:27  h  47.40  2. 

ML522126.IMG  5/22/95  10:50:31  v  47.40  2, 

ML522127.IMG  5/22/95  10:50:58  u  44.70  2. 

ML522128.IMG  5/22/95  10:51:11  h  47.10  2. 

ML522129.IMG  5/22/95  10:51:16  v  47.10  2. 

ML522130.IMG  5/22/95  10:51:27  u  44.90  2. 

ML522131.IMG  5/22/95  10:53:00  h  47.10  2. 

ML522132.IMG  5/22/95  10:53:05  v  47.10  2. 

ML522133.IMG  5/22/95  10:53:14  u  43.70  2. 

ML522134.IMG  5/22/95  10:53:27  h  46.90  2. 

ML522135.IMG  5/22/95  10:53:32  v  46.90  2. 

ML522136.IMG  5/22/95  10:53:41  u  44.40  2. 

ML522137.IMG  5/22/95  10:54:51  h  47.10  2, 

ML522138.IMG  5/22/95  10:54:58  v  47.10  2. 

ML522139.IMG  5/22/95  10:55:07  u  44.20  2. 

ML522140.IMG  5/22/95  10:55:21  h  47.10  2. 

ML522141.IMG  5/22/95  10:55:25  v  47.10  2. 

ML522142.IMG  5/22/95  10:55:33  u  44.70  2. 

ML522143.IMG  5/22/95  11:05:25  h  47.10  2. 

ML522144.IMG  5/22/95  11:05:30  v  47.10  2. 

ML522145.IMG  5/22/95  11:05:42  u  44.40  2. 

ML522146.IMG  5/22/95  11:05:55  h  47.60  2. 

ML522147.IMG  5/22/95  11:06:00  v  47.60  2. 

ML522148.IMG  5/22/95  11:06:10  u  44.40  2. 

ML522149.IMG  5/22/95  11:08:39  h  47.40  2. 

ML522150.IMG  5/22/95  11:08:45  v  47.40  2. 

ML522151.IMG  5/22/95  11:08:54  u  44.20  2. 

ML522152.IMG  5/22/95  11:09:10  h  47.40  2. 

ML522153.IMG  5/22/95  11:09:15  v  47.40  2. 

ML522154.IMG  5/22/95  11:09:25  u  44.40  2. 

ML522155.IMG  5/22/95  11:11:18  h  47.90  2. 

ML522156.IMG  5/22/95  11:11:23  v  47.90  2. 

ML522157.IMG  5/22/95  11:11:34  u  44.90  2. 

ML522158.IMG  5/22/95  11:11:52  h  47.60  2. 


119 


4. 

0. 

0. 

4. 

0. 

0. 

4. 

703. 

3311 

4. 

0. 

0. 

4. 

0. 

0. 

4. 

894. 

2818 

4. 

0. 

0. 

4. 

0. 

0. 

4. 

156. 

5117 

4. 

0. 

0. 

4. 

0. 

0. 

4. 

0. 

0. 

4. 

4. 

6647 

4. 

0. 

0. 

4. 

0. 

0. 

4. 

5. 

6955 

4. 

0. 

0. 

4. 

0. 

0. 

4. 

690. 

3843 

4. 

0. 

0. 

4. 

0. 

0. 

4. 

1143. 

2748 

4. 

0. 

0. 

4. 

0. 

0. 

4. 

434. 

4618 

4. 

0. 

0. 

4. 

0. 

0. 

4. 

436. 

4469 

4. 

0. 

0. 

4. 

3. 

0. 

4. 

365. 

5244 

4. 

4. 

0. 

4. 

0. 

0. 

4. 

1. 

6435 

4. 

0. 

0. 

4. 

0. 

0. 

4. 

1699. 

1756 

4. 

0. 

0. 

4. 

0. 

0. 

4. 

913. 

3436 

4. 

0. 

0. 

4. 

0. 

0. 

4. 

1175. 

3080 

4. 

0. 

0. 

4. 

0. 

0. 

4. 

465. 

4504 

5. 

0. 

0. 

5. 

905. 

0. 

5. 

1680. 

4428 

5. 

0. 

1. 

5. 

0. 

0. 

5. 

1569. 

5269 

5. 

0. 

0. 

5. 

1141. 

0. 

5. 

1784. 

4135 

5. 

32. 

0. 

5. 

1116. 

0. 

5. 

2206. 

2276 

5. 

0. 

3. 

5. 

0. 

0. 

5. 

1362. 

6366 

5. 

0. 

1171 

ML522159.IMG  5/22/95  11:11:59  v  47.60  2 

ML522160.IMG  5/22/95  11:12:09  u  44.40  2 

ML522161.IMG  5/22/95  11:19:45  u  44.40  5 

ML522162.IMG  5/22/95  11:22:19  h  47.90  2 

ML522163.IMG  5/22/95  11:22:24  v  47.90  2 

ML522164.IMG  5/22/95  11:22:54  u  44.20  2. 

ML522165.IMG  5/22/95  11:23:17  h  48.10  2 

ML522166.IMG  5/22/95  11:23:25  v  48.10  2, 

ML522167.IMG  5/22/95  11:23:42  u  44.20  2, 

ML522168.IMG  5/22/95  11:25:56  h  47.60  2. 

ML522169.IMG  5/22/95  11:26:04  v  47.60  2, 

ML522170.IMG  5/22/95  11:26:19  u  44.40  2. 

ML522171.IMG  5/22/95  11:26:34  h  48.10  2. 

ML522172.IMG  5/22/95  11:26:41  v  48.10  2. 

ML522173.IMG  5/22/95  11:27:16  u  44.40  2. 

ML522174.IMG  5/22/95  11:28:27  h  47.60  2. 

ML522175.IMG  5/22/95  11:28:32  v  47.60  2. 

ML522176.IMG  5/22/95  11:28:45  u  44.20  2. 

ML522177.IMG  5/22/95  11:28:55  h  47.60  2. 

ML522178.IMG  5/22/95  11:29:06  v  47.60  2. 

ML522179.IMG  5/22/95  11:29:23  u  44.40  2. 

ML522180.IMG  5/22/95  11:38:45  h  47.60  2. 

ML522181.IMG  5/22/95  11:38:50  v  47.60  2. 

ML522182.IMG  5/22/95  11:39:32  u  44.20  2. 

ML522183.IMG  5/22/95  11:40:07  h  47.90  2. 

ML522184.IMG  5/22/95  11:40:13  v  47.90  2. 

ML522185.IMG  5/22/95  11:40:40  u  44.20  2. 
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ML523022.IMG  5/23/95   8:47:04  u  43.50  2. 
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APPENDIX  F.   IDL  IMAGE  PROCESSING  CODES 

The  following  programs  are  used  to  determine  TBRD 
automatically   from  a  file  which  contains  all  the  specific 
information  about  a  target.   The  desired  information  is 
shown  in  the  file  compilation  of  image  names  in  Appendix  E. 

Also  the  IDL  code  to  subtract  the  two  split-images  is 
contained  in  this  appendix.   The  file  headers  denote  the  use 
of  the  code  and  the  code  itself  is  commented  to  provide  its 
own  detailed  explanation. 
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;fliral2 .pro  processes  a  single  ship  for  contrast  using  sreadpol  to  get 

filenames  from 

; experiment  flirl2 

;Pete  version  that  uses  new  contrast  with  no  denominator 

; taken  from  flira8  which  selects  files  directly 

; 2-17-95  revised  xO,ix  for  xO  always  to  be  peak,  changed  shipbox  to  find 

peak  if  the 

;new  peak  selected  with  the  cursor  is  eliminated  by  the  threshold 

;xO,yO         peak  value  in  contiguous  area 

;ix,iy,nx,ny   box  coordinates 

; ship  temperature  ship 

; shipb         raw  data  ship 

;yc  horizon  line  point 

;hz  iy-yc    location  of  divisor  in  ship  box 

.*********************************************************************** 

******************* 

; routines  readflir, shipbox, thresh, remake, fileonly, sreadpol  used 

;readflir  reads  file  data  from  flir  format 

; horizon  allows  choice  of  horizon,  formerly  profiles 

; shipbox  finds  contiguous  ship  and  makes  ship  box 

;threshld  finds  and  returns  edge  threshold  thresh  from  yy  and  t  as  well 

as  peak  xO,yO 

; remake  uses  threshold  threshO  to  make  ship  and  shipb  from  t  and  yy 

; fileonly  separates  filename  from  file  and  path  for  output 

; sreadpol  reads  polarization  file  and  separates  h  or  v  in  variable  pol 

!order=0  ; get  rid  of  old  program 

inverision 

iscale=2  ; display  2x  size 

on_ioerror, problem  ;if  error,  quit  program 

************************************************************************ 

******************** 

; select  the  files  with  wild  cards 

pathl=' ' 

file=pathl 

if (n_elements (path)  ne  0)  then  print, ' oldpath=  ',path 

read, 'input  path  for  files  or  cr  for  oldpath' ,pathl 

if(pathl  ne  ' ' )  then  path=pathl+* \ '   else  path=* . * 

if (strmid (path, strlen (path) -1, 1)  ne  ' \')  then  path=path+' \ '         ; f ix 

oldpath  error 

read, 'input  filename  of  polarization  information' , file 

sreadpol, file, files, pol  ;get  files  from  flirllp 

made  separately 

files=path+files  ;put  files  in  standard 

findfile  format 

filelist, files, path  ;program  to  list 

filenames  that  are  different 

; f or  i=0,num-l  do  print, files (i) , i  ;print  the  filenames  with 

a  number 

i='  ' 

RESTART:  read, ' Input  start  and  stop  number  for  processing  CR  to  stop',i 

if  strlen  (i)  eq  0  then  goto,  finish 

istart=f ix  (i) 

istop=istart 

read, 'Stop  number    ' , istop 


126 


C5N=DBLARR( ( ISTOP-ISTART+1 ) ,3)  ;The  first  z  is  full  box,  second,  above 

horiz,  third  below 

read, 'input  first  guess  difference  t  level,  eg  1  to  include  all  of 

ship' , tdif 
.*********************************************************************** 

******************** 

close, 2 

fileout='  ' 

read, ' input  output  filename  for  contrast  results ', fileout 

openw, 2 , fileout 

printf , 2, systime ( ) , '   date  of  analysis' 

printf , 2, istop-istart+1, '    number  of  files' 
.*********************************************************************** 

r 
******************** 

auto='      * 

read, 'input  CR  for  semi-automatic  proccessing' ,  auto 

tdif save=tdif  ; save  tdif  so  that  changes  will 

not  affect  batch 

printf , 2, 'Temperature  difference  =  ',tdif  ; save  in  a  file 

; tdif ,nx,ny, yc, iy 

printf, 2, ' filel, tdif , nx,ny, yc (horizon) , iy (ship  bottom) , Polarization' 

for  i=istart,  istop  do  begin        ;ii  =  first  of  pair 

tdif=tdif save  ; restore  the  temperature  to  the 

original  batch  value 

f=files(i) 

f ile=f  ; save  the  name 

readflir, f, yy, t, bar, tt, head         ; read  flir  data,  convert  to 

temperature,  add  bar 

-  *****************************cjisDlav  the  raw 
(-j  g  •(-  ^  ****************************************  * 

loadct,0  ; set  color  table  0 

iscale=2  ;make  pictures  twice  their  normal 

size 

; remove  the  path  length  from  the  file  string  and  label  the  tv 
image  with  the  filename 

pathlen=strlen (path) +1 

filelen=strlen ( file) 

labelen=filelen-pathlen 

label=strmid ( file,pathlen, labelen) 
window, 0, xsize=140*iscale, ysize=140*iscale+50, xpos=0, ypos=0, title=' 0 
*+label 

tvscl, rebin (tt, 140*iscale, 165*iscale)   /display  the  image 
.  ********************************************£^1(-j  of  temperature 

Q^g-Q]_£y*********  ************ 

window, 2,xpos=0, ypos=165*iscale+35, xsize=140*iscale, ys i ze=l 40 *i scale 

print, ' choose  horizon  with  left  button  or  reject  picture  with  right' 

tvscl, rebin (yy, 140*iscale, 140*iscale)  ;display  the 

binary  image  raw 

HORIZON, smooth (rebin (yy, 140*iscale, 140*iscale) , 5) ,xc, yc     ; read  cursor 

position 

if(!err  eq  4)  then  goto,  skiploop  ;skip 

processing  of  file 

xc=xc/iscale 

yc=yc/iscale  ; scale  to  140 

size  not  picture  size 

print , ' xc , yc=  ' , xc , yc 

;  ****************************t}ireshold  with   temperature 

cji  ffg^gi^ce*************************** 
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threshld, t, yy, thresO, xO, yO  ; define 

thres  and  find  xO,yO  peak 

;RESHIP  RESTORES  THE  VARIABLES  MODIFED  BY  SHIPBOX,  NOT  THRESH 

errorct=0  ; don't  allow  infinite 

loopes  of  choices 

***************************************  *RESHI  p************************** 

********************* 

reship: 

remake, ship, shipb, t, yy, thresO  ; remake  the 

thresholded  ship  and  shipb  using  thresO 

tvscl, rebin (shipb, 140*iscale, 140*iscale)  ;display  the  binary 

image  raw 

tvcrs,xO*iscale, yO*iscale  ; cursor  on  hot  spot 

shipbox, shipb, ship,xO, yO, ix,iy,nx,ny, iscale        ; ix, iy,nx,ny  define 

ship  box  contiguous 

;   WARNING  ALL  VARIABLES  EXCEPT  ISCALE  ARE  CHANGED  IN  THIS 
CALL 

boxdraw, ix*iscale,iy*iscale,nx*iscale,ny*iscale     ; show  where  the  box 
is  in  *iscale 

tvcrs,xO*iscale, yO*iscale, /device  ;put  the  cursor  on 

the  peak  location 
************************************************************************ 

********************* 

;***above  makes  box  with  minimum 

•(-J-^ggJ-jO^,-^*************************************************** 

! ERR=4 

; if (auto  ne  ' *  or  errorct  gt  0)  then  begin  /suspend  auto  only 

for  errors 

print, 'right  button,  use  auto  position,  left,  new  cursor  position' 
cursor, ixnew, iynew, /device, /down  ;when  pushed, 

choose 

print, ' !ERR=  ', ! ERR 

if(!ERR  eq  4)  then  print, 'auto  selected' 
if(!ERR  eq  1)  then  begin 
x0=ixnew/ iscale 
y0=iynew/ iscale 

print, 'new  peak  at  ',x0,y0  ;NEW  peak 

at  x0,y0 

remake, ship, shipb, t, yy, thresO  ;if  not 

remade,  the  shipbox  will  have  destroyed  shipb 

; WARNING  SHIPBOX  DESTROYS  SHIPB  AND  SHIP 

shipbox, shipb, ship,x0, y0,ix, iy,nx,ny, iscale         ;ix,iy,nx,ny 
define  ship  box  contiguous 

boxdraw, ix*iscale,iy*iscale,nx*iscale,ny*iscale     ; show  where 
the  box  is  in  *iscale 

endif 
;endif  ; semi  auto  is 

two  button  pushes  now 

.*****************IF  THE  H0T  SPOT  IS  FOUND  DO  LOCAL  THRESHOLD  IN  THE 

BOX** ************* 

remake, ship, shipb, t, yy, thresO+tdif  ;  remake  with  full 

threshold 

print, ' redisplay  with  full  threshold' 

shipbox, shipb, ship, xO, yO, ix, iy,nx,ny, iscale  ;x0,y0,nx,ny 

define  ship  box  contiguous 

tvscl, rebin (shipb, 140*iscale, 140*iscale) 

boxdraw, ix*iscale, iy*iscale,nx*iscale, ny*iscale 
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tvcrs,ix,iy  ;place  the  cursor  on 

the  max  for  a  check 

print, ' cursor  on  maximum1 

. ************************* fi^   a  new  shipbox  based  on  threshold  and 

tdif*  ****************** 


;if  data  is  all  zero, 
; avoid  an  infinite 


if  (total (shipb)  eq  0)  then  begin 
threshold  did  not  work 

errorct  =  errorct+1 
loop 

print, 'no  ship  was  found.  Data  is  not  usable  with  this  tdif',tdif 

if  (errorct  gt  2)  then  begin 
print,'  skipping  ' , f ile 
goto,  skiploop 

endif 

read, 'input  new  tdif',tdif 
;       auto='  '  ;  suspend  auto 

permanently  on  first  error 

goto,  reship 

endif 

tdiff=' ' 

print, 'tdif=  *,tdif 

tdiff=' ' 

if (auto  ne  ' '  or  errorct  gt  0)  then  $ 

read, 'input  new  tdif  or  CR  for  old  or  negative  to  skip',tdiff 

if (strlen (tdif f }  eq  0)  then  tdiff=tdif  else  begin 

tdif=float(tdiff) 

if (tdif  It  0)  then  goto,  skiploop 

goto, reship 
endelse 


-*•********************************************************************* 

r 
******************** 

; choices  section  checks  to  see  which  picture  is  larger  and  uses  the 

appropriate  tdif 

; to  remake  and  replot  the  temperatures  reference  windows 

; ship  temperature  ship  in  box    &  sv  extension 

;box    temperature  box  not  ship 

; shipb  binary  ship  in  box 
;boxb  binary  box 


ship=(t  gt  thresO  +  tdif ) *t 
the  OLD  threshold 

shipb=(t  gt  thresO  +  tdif)*yy 


;be  sure  NOR  to  use 


ship=ship (ix: ix+nx, iy: iy+ny) 
shipb=shipb (ix: ix+nx, iy : iy+ny) 
box=t (ix: ix+nx, iy : iy+ny) 

box  from  original  data 

boxb=yy ( ix : ix+nx, iy : iy+ny) 
box=box* (ship  eq  0) 

(temperature)  in  box 


; subset 
/everything  in  the 

;not  ship 
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boxb=boxb* (shipb  eq  0) 
zscale=140/nx 

window, 5,xpos=140*iscale,ypos=0,xsize=nx*zscale, ysize=ny*zscale 
tvscl, congrid (ship, nx*zscale,ny*z scale) 
-*********************************************************************** 

******************** 

;make  contrasts   ell  =  first  file,  shipfull,  cl2=above  horizon,  cl3=up 

to  horizon 

;C3N   full  box,  above  horizon,  below  horizon 

; yc  horizon  line  point 

; iy  ship  bottom 

;hz  yc-iy    location  of  divisor  in  ship  box 

.*********************************************************************** 

******************* 

hz=yc-iy 

;if (hz  It  0)  then  horizon  is  below  ship  and  ship  is  flying 

;if(hz  ge  iy+ny)  then  horizon  is  completely  above  ship 

c5n (i-istart, 0) =pcontra ( shipb, boxb) 

; f rom  yc (horizon)  -iy  to  ma 

if (hz  le  0  or  hz  ge  ny)  then  begin  ;hz  0  means  c5n(i,2) 

is  ill  defined  as  ship  is  on  horizon 

print,'hz,ny  ',hz,ny, '  error'  ;hz  gt  ny  means 

c5n(i,l)  is  below  horizon 

c5n (i-istart, 1) =c5n (i-istart, 0)  ;hz  It  0  means  ship 

is  flying  above  horizon 

c5n (i-istart, 2) =c5n (i-istart,  0) 
endif  else  begin 

c5n (i-istart, 1) =pcontra (shipb (*,hz : * ) ,boxb (*,hz: *) ) 

c5n (i-istart, 2) =pcontra (shipb (*, 0 :hz) ,boxb (*, 0 :hz) ) 

; duplicate  the  first  if  the  value  hz  is  zero 

endelse 

print, fileonly (files  (I ), path) , '  ' , tdif ,nx,ny, yc, iy, '  ',pol(i) 
print, C5N( (I-ISTART) ,0) , f ormat="  (  (F7.3,2x) ,  'Full  ship' ) " 
print, C5N( (I-ISTART) ,1) , f ormat=" ( (F7.3,2x) , 'Above  Horizon' ) " 
print, C5N( (I-ISTART) ,2) , f ormat=" ( (F7.3,2x) , 'Below  Horizon' ) " 

print f, 2, fileonly (files (I) ,path) , '  ' , tdif ,nx,ny, yc, iy, '  ' ,pol (i) 
printf,2,C5N( (I-ISTART) ,0) , f ormat=" ( (F7.3,2x) , 'Full  ship' ) " 
printf,2,C5N( (I-ISTART) ,1)  ,  f ormat=" ( (F7.3,2x) , 'Above  Horizon' ) " 
printf,2,C5N( (I-ISTART) ,2) , f ormat=" ( (F7.3,2x) , 'Below  Horizon' ) " 


skiploop: 

endfor 
.*********************************************************************** 

r 
******************** 

goto, finish 

problem:  print, 'I  had  a  problem  and  quit' 

finish:  print, 'program  finished' 

close, 2 

end 

; program  shipszl2  reads  output  of  fliral2  (Pete  Heisey)  which  is: 
/filename,  tdif ,ix, iy,nx,ny,H  or  V 
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;and  calculates  TBRD  with  ship  size  and  plots  it 

close, 1  ; shut  up  in  case  of  an  error 

tempp= ' 

file=' 

read, ' input  filename  from  analysis  output1, file 

on_ioerror, stoperr 

openr, 1, file 

nf iles=0 

iy=0 

yc=iy 

nx=iy 

ny=iy 

i=0 

t-1.0 

HV=*  '  ;file  count 

while  not  eof(l)  do  begin  ; read  to  end  of  file 

which  is  not  more  than  nfiles 

;loop: readf , 1, tempp  &  print, tempp  ; read  string  and  print 

it 

loop: readf, 1, tempp 

if (nfiles  eq  0)  then  begin  ;if  number  of  files  not 

fount  yet,  look  for  it 

if (strpos (tempp, 'number  of  files')  gt  0)  then  begin 

reads, tempp, nfiles  &  print, tempp, nfiles      ; get  the  number  of  files 

rdata=fltarr (5, nfiles)  ;make  array  of  number  of 

files,  may  not  use  it  all 

cdata=fltarr (nfiles, 3)  ;array  of  contrast  data 

;    namearr=strarr (12, nfiles)  ;  array  of  filenames 

namearr=strarr (nfiles) 

pol=strarr (nfiles) 

print, ' finished  making  arrays,  Nfiles=  '  , nfiles  ;let  me 

know  arrays  were  made 

endif 
endif 
if (strpos (tempp, '. IMG' )  It  0)  then  goto,  loop  else  $ 

namearr (i) = (strmid (tempp, 0, strpos (tempp, ' G' ) +1) )  ; search  for  file  end 
IMG 

;print,  (strmid (tempp, strpos (tempp, ' G' ) +1,  strlen (tempp) ) ) 

y=strcompress (strmid (tempp, strpos (tempp, ' G' ) +1, strlen (tempp) ) )      ; get 

rid  of  white  spaces 

; reads, y, t, ix,iy,nx,ny, HV  &  print, t, ix, iy, nx, ny,HV    ; convert  data  and 

print  it 

;printf , 2, tdif ,nx,ny, yc, iy  maee  this 

reads, y,t,nx,ny,yc,iy,HV 

rdata (0,i)=t 

rdata (1, i) =nx 

rdata (2,i) =ny 

rdata (3, i) =yc 

rdata (4, i) =iy 

pol (i)=HV 

y=1.0 

readf, l,y  &  cdata(i,0)  =  y 

readf, l,y  &  cdata(i,l)  =  y 

readf, l,y  &  cdata(i,2)  =  y 

print, ' incrementing  loop  counter' 
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i-i+1 

variable 

if (rdata (4,i-l)  le  1)  then  i=i-l 

the  data 

endwhile 

goto,  stopp 

stoperr :print, ' io  error' 

goto, loop 

stopp: close,  1 

if(i  It  0)  then  i=0 

print, ' finished1 , i 

rdata=rdata (*,0:i-l) 

cdata=cdata (0 : i-1, *) 

namearr=namearr (0 : i-1) 

print,' file=  ' , f ile 

S PLOT 9e, RDATA, CDATA, FILE, pol 

end 


; increment  the  stored 
;if  height  is  0,  skip 


/program  contrast. pro  which  calculates  TBRD  for  each  ship  image 
;anayzed.   it  is  a  subroutine  of  shipszl2 

function  contrast, ship, box 
contrast= (shipav (ship) -ship (box) ) 
return, contrast 
end 


;program  splot  9  is  subroutine  of  shipszl2  which  takes  the  ship 
conatrast  data  and  plots  it  vs .  ship  range  determined  by  size 

PRO  SPLOT9, RDATA, CDATA, file, pol 

ptitle=file 

;cdata(*,0)     contrast  for  whole  ship 

;cdata(*,l)     contrast  for  above  horizon 

;cdata(*,2)     contrast  for  below  horizon 

; rdata (2,*)     ship  size  from  height 

hpol=where (strtrim(pol, 2)  eq  'h')       ; array  of  number  for  h 

polarization 

vpol=where (strtrim(pol,  2)  eq  *v* )       ;array  for  v  polarization 

help, hpol , vpol 

if(n  elements (hpol)  le  1)  then  hpol=where (strtrim(pol, 2)  ne  'h') 
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if (n_elements (vpol)  le  1)  then  vpol=where (strtrim(pol, 2)  ne  'v') 

LOADCT,39 

; read, ' Input  title ' ,ptitle 

;WINDOW,0,XSIZE=1024, YSIZE=768 

;plot, l./rdata(2, *) , cdata (*,2) , CHARSIZE=2,psym=2, color=0,background=255, 

charthick=2, title=' ! 17'+ptitle  $ 

; , xtitle= ' Ship  Range ', ytitle=' Target-Background  Radiance 

Difference ' , /nodata 

; oplot, 1 . /rdata (2, * ) , cdata (*,  2 ) , psym=2, color=0 

; oplot, 1 . /rdata (2, * ) , cdata (*, 0) ,psym=l, color=0 

; oplot, 1 . /rdata (2, *) , cdata (*, 1) ,psym=4, color=60 

;xyouts,  .65,  .87, charsize=3, charthick=2, ' ! 17Above 

Horizon' , color=60, /normal 

;xyouts, .65, .82, charsize=3, charthick=2, ' ! 17Full  Ship' , color=0, /normal 

;xyouts, .65, . 77, charsize=3, charthick=2, ' ! 17Below 

Horizon ' ,  color=60, /normal 

;WINDOW, 1,XSIZE=1024,YSIZE=768 

;plot, 1. /rdata (2,*) , cdata (*, 0) , CHARSIZE=2,psym=2, color=0,background=255, 

charthick=2,title=' !17'+ptitle  $ 

; ,xtitle='Ship  Range ', ytitle=' Target-Background  Radiance  Difference1 

; oplot, 1. /rdata (2, * ) , cdata (*, 1) ,psym=4, color=0 

; ; oplot,  1 . /rdata (2, * ) , cdata (*,  2 ) ,psym=6, color=0 

WINDOW,2,XSIZE=1024,YSIZE=768 

plot, 1. /rdata (2, *) , cdata (*,2) , CHARSIZE=2,psym=2, color=0,background=255, c 

harthick=2,title=' ! 17 ' +ptitle, /nodata  $ 

,xtitle=' Ship  Range  from  Ship  Size* , ytitle=' Target-Background  Radiance 

Difference' 

oplot , 1 . / rdata ( 2 , hpol ) , cdata ( hpol , 2 ) , psym=2 , color=0 

oplot, 1. /rdata (2, vpol) , cdata (vpol, 2) , psym=6, color=0 

oplot, 1 . /rdata (2, hpol) , cdata (hpol, 1) ,psym=l, color=0 

oplot, 1 . /rdata (2, vpol) , cdata (vpol, 1) ,psym=4,  color=0 

xyouts,  .65,  . 77,charsize=l,charthick=2, ' ! 17  *+  =  H   ! 9B  ! 9V  ! 17= 

V , color=0, /normal 

xyouts, . 65,  . 87, charsize=l, charthick=2, '! 17Above  Horizon   * 

! 9B' ,color=0, /normal 

xyouts, . 65, . 82, charsize=l, charthick=2, '! 17Below  Horizon  + 

! 9V ,color=0, /normal 

end 


;  this  program  subtracts  the  arrays  of  the  two  split  images  from  the 

splitimage  polar imeter. 

;program  splitcom  compares  split  image  in  yy  array 

getfiles, files 

; get  array  of  filenames  in  files 
loop: 

;do  it  again 
read,  'input  number  of  file  to  get',ifile 
if(ifile  It  0)  then  goto, problem 
f=files(ifile) 

readflir, f, yy, t, bar, tt, head  ; get 

file  and  make  temperature  array 
! order=0 
on_error, 1 
ybot=yy(0:*,70:*) 
ytop=yy (0:*,0: 69) 
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read, 'input  y  to  invert ', yes 

if (yes  eq  ' y')  then  begin 

ybot=ybot (*, 69-indgen (70) )  ; Invert  the  image 

ytop=ytop (*, 69-indgen (70) ) 

endif 

hotspot, ytop,xhott, yhott  ; FIND  THE 

HOTS POT  IN  THE  IMAGE 

hotspot, ybot, xhotb, yhotb  ; subroutine 

window, l,xsize=140*3,ysize=140*3 

tvscl, rebin (ybot, 140*3,70*3) 

tvscl, rebin (ytop, 140*3, 70*3) , 0, 70*3        ;display  both  images 

xdif=xhott-xhotb 

/calculate  displacement  to  overlay  hotspots 

ydif=yhott-yhotb 

print,  xdif,  ydif,  '  xdif  ydif 

ytop=float (ytop  )  /make 

floating  so  subtraction  will  work 

ybot=float (ybot) 

displace, ybot, fix (xdif ), fix (ydif )  ;  subroutine  which 

; overlaps  images 
window, 3, xsize=140*4, ysize=70*4 

absoldif=abs (ytop- ybot)  ;make  the 

values  of  the  difference  positive 

tvscl, rebin (absoldif, 140*4,70*4)  /display  difference 

goto, loop 
problemrprint, 'exiting' 

/error  and  normal  exit  for  -1  filename 
end 


/pro  di spl ace, yval, xdif, ydif 

ysize=size (yval) 

print, ysize 

if  (ydif  ge  0)  then 

yval (0:*,abs (ydif) : ysize (2) -1) =yval (0 : *, 0 : ysize (2) -1-abs (ydif) )  $ 

else  yval (0 :*, 0: ysize (2) -1-abs (ydif) ) =yval (0 : *, abs (ydif ) : ysize (2) -1) 
end 

pro  hotspot , yval , xhot , yhot 

ysize=intarr (5) 

ysize=size (yval) 

print, max (yval, ymax)   &  print, ymax 

print,  ymax- (ymax/ysize (1) ) *ysize (1)  , ymax/ysize (1) 

xhot=  ymax- (ymax/ysize (1) ) *ysize (1) 

yhot=ymax/ysize (1) 

end 
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APPENDIX  G. 
SPLIT  IMAGE  POLARIMETER  IMAGES,  3  AUGUST  1995 


The  following  appendix  contains  the  images  recorded  of 
the  USS  John  McCain  (Arleigh  Burke  class  DDG)  on  3  August, 
1995  with  the  split-field  polarimeter.   This  was  a 
preliminary  experiment  conducted  to  test  the  performance  of 
the  optical  system  in  the  far  field  and  to  validate  field 
performance  of  all  system  components.  Any  Data  on  the  image 
file  not  contained  in  the  Experimental  notes  is  stored  in 
the  image  file  header. 
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Experimental   Notes 


Image  File 
Name 

Targe 

t 

Thermal 
Range 

Remarks 

m0803001 

BB 

setup  test-BB=black  body  source 

m0803002 

BB 

infinite  focus,  BB  @20' 

m0803003 

BB 

ii 

m0803004 

BB 

Black  body  temperature  36degrees 

m0803005 

McCain 

2 

cal  table  4003  used 

m0803006 

McCain 

5 

m0803007 

McCain 

2 

m0803008 

McCain 

2 

m0803009 

McCain 

2 

m0803010 

McCain 

5 

m0803011 

McCain 

2 

m0803012 

McCain 

5 

m0803013 

BB 

2 

cal  tab.  4004,  BB  36  degrees 
OC=-4  4 

m0803014 

BB 

5 

m0803015 

McCain 

5 

m0803016 

McCain 

2 

m0803017 

McCain 

2 

m0803018 

McCain 

2 

m0803019 

FOV 

2 

empty  FOV  with  aft  of  ship  on 
far  right  side 

m0803020 

McCain 

5 
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m0803021 

McCain 

10 

m0803022 

McCain 

5 

m0803023 

BB 

10 

focus  set  on  inf. 

for  BB 
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